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CHAPTER I 


SUMMARY 


1.0 INTRODUCTION 


The two principal objectives of this report are (1) to 
describe the NMSU/PSL Radiometer Antenna Calibration Facility 
(RACF) and (2) to summarize the antenna and radome loss measure- 
ments made on the Passive Microwave Imaging System (PMIS) and 
the Multifrequency Microwave Radiometer (MFMR) during January 
and February of 1975, This chapter summarizes the major 
features of the facility, points out highlights of the PMIS/MFMR 
loss measurements, and makes several recommendations for future 
measurements . 

1.1 NMSU/PSL Radiometer Antenna Calibration Facility 


The physical facility used at PSL for the loss calibration of 
antennas and radomes consists of a large reflecting bucket and 
an electronic equipment building, both located at A-Mountain, 
approximately 5 km east of the New Mexico State University. The 
bucket is at an altitude of 1.46 km above MSL and is an aluminum 
foil covered wooden truncated inverted pyramid with dimensions as 
shown in Fig. 2-7. The purpose of the bucket is to block thermal 
emission from surrounding terrain and from near-horizon atmospheric 
sources. Thus, the thermal emission incident on the antenna is 
ideally isothermal and equal to the zenith sky temperature at 
the frequency of interest. 

In addition to the physical plant, an extensive data 
reduction program has been written for PMIS and MFMR, with 
separate software components tailored to the PCM output format 
used in the NASA P-3A earth resources aircraft. This is discus- 
sed further in Chapter III and the Appendix. 


2 


The bucket is similar to an enclosure at Table Mountain, 
California used previously by NASA. However, the timely 
availability of PSL support personnel and a real-time data 
reduction capability make the RACF more useful as a radiometer 
test bed. 

1.2 Highlights of Test Results 

It is shown conclusively in Chapter VIx that the bucket 
technique provides sufficient accuracy and repeatability for 
the loss values to be used in an aircraft radiometric measurement 
mission when the radiometer aircraft mockup exhibits sufficient 
fidelity to the actual aircraft situation. An extensive error 
analysis for PMIS and MFMR (Secs. 5.5 and 6.7) establishes the 
following: 

1. For PMIS (10.69 GHz) the total uncertainty (sum of 
random plus systematic errors) in the antenna loss is + 0.022 

and + 0.039 for the vertical and horizontal channels respectively. 
The total uncertainty in the radome loss is + 0.014 and + 0.039 
for the vertical and horizontal channels. This would correspond 
to a total uncertainty in the measured brightness temperature 
(in flight) of +5.2 K and +5.7 K for the vertical and horizontal 
channels respectively. 

2. For MFMR, L-Band (1.4135 GHz) , the bulkhead and associated 
absorber were not sufficiently good replicas of the actual aircraft 
situation, with the result that the mutual coupling between the 
AIL array and the radome was so strong that the measured values 

of the radome loss were virtually meaningless. However, as 
explained in Sec. 6.6, a constant value of L R = 1.09 is a reasonable 
choice. The antenna loss values were typically « 1.36 (es- 
sentially independent of pitch angle) when the PSL-furnished 
X-Band absorber was in place on the bulkhead. However, this 
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material is virtually transparent at L-Band and thus does not 
simulate the actual aircraft situation where an L-Band absorbing 
material is used. Thus, in the flight situation, it may be 
expected that L_ > 1.36 for pitch angles near 0° or 180°, since 
the relatively hot absorber will contribute noise power through 
the antenna sidelobes. In the absence of accurate near-field 
in situ patterns for the AIL array, it is difficult to see how 
the effect of the absorber can be estimated, short of new 
improved measurements. 

3. For MFMR, K , K, K -Bands (18.0, 22.05 and 37.0 GHz), 

U cL 

the antenna loss is less than the radome loss and both losses 
depend much more on the jpitch and roll angles than was the case 
at L-Band. The antenna loss is slightly higher in Channel 1 
than in Channel 2 and generally rises as the pitch angle approaches 
180° when the horns are nearest the bulkhead (c.f. Figs. 6-2, 6-3, 
and 6-4). However, as explained in Sec. 6.5), the exact values 
of L^ and the slope dL^/d0 Q depend critically on the fidelity of 
the mockup used, the exact shape of the absorber used on the 
bulkhead, etc. 

The radome loss at K^, K, and K a ~Bands also depends critically 

on both pitch and roll angles, as shown in Figs. 6-5 - 6-10. The 

jagged nature of these loss graphs is not due to errors in the 

measurement (the repeatability was much better than this) , but 

seems to stem from resonant scattering between the horns and 

radome material or radome resonant thickness effects (especially 

at K -Band) . 
a 


Thus, it may be stated that the precision of measurement was 
good, but the accuracy was poor, since the antenna system furnished 
was an insufficiently good replica of the antenna system used. 

The MFMR error budget is discussed at considerable length in 
Sec . 6.7. 
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1.3 Recommendations 


As a result of the measurement program, several major 
recommendations can be made: 

1. Future tests of the MFMR should be made with a much 
better replica of the aircraft bulkhead, including the same 
absorber that would be used in flight operations, along with any 
other significant structural features such as weather radars, etc. 

2. Measurements of the input vs. pitch angle at L, K u , K, 

and K -Bands should be made, with and without the radome. 
a 

3. A study of the MFMR pitch and roll positioning repeat- 
ability should be made. 

4. If the K , K, and K -Band channels are to be used for 

u f a 

skyward viewing, the horns should be located at a roll angle of 
180°, rather than 0°. This would obviate the critical dependence 
on the exact bulkhead geometry. 

5. Greater attention should be paid to reducing PGM noise 
and long-term instability in the PMIS receivers, particularly the 
horizontal channel. 


5 


CHAPTER II 


ANTENNA LOSS MEASUREMENT TECHNIQUES 
2 . 0 INTRODUCTION 


The temperature calibration of a radiometer system normally 
proceeds in two parts: (1) the calibration of the connecting 

waveguide, front end, etc., by use of known reference load tempera- 
tures, and (2) the calibration of the antenna and radome loss. The 
purpose of the NMSU/PSL Radiometer Antenna Calibration Facility 
is to measure the antenna and radome loss by using a calibrated 
radiometer receiver and a known source brightness temperature. 

This facility differs by comparison to many other past cali- 
bration efforts, in that relatively large aircraft and satellite- 
borne radiometer systems, including automatic positioners and 
radomes, may be accommodated. Also, the bucket technique 
employed is easily adaptable to multi- frequency use since radio- 
sonde support is locally available on a timely and convenient 
basis. The genesis of this effort was the requirement to 
periodically calibrate the NASA JSC Passive Microwave Imaging 
System (PMIS) at 10.69 GHz and the Multifrequency Microwave 
Radiometer (MFMR) at 1.4, 18, 22.05, and 37 GHz. These systems 
are complicated by a variety of requirements for beam scanning, 
polarization switching, multiplexed data handling and subsequent 
data reduction. 

In response to the need for near real-time data reduction, 

PSL has developed an extensive set of software packages to 
handle the PCM encoded radiometer and housekeeping data, both for 
PMIS and MFMR. These make use of a PCM decommutation device and 
an IBM System 7 computer which transfers the stripped data from 
the magnetic tape to an IBM System 370 Computer for processing. 

This will be described in more detail in Chapter 3. 
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2 . 1 Definitions 


The antenna loss of a receiving antenna is defined as 

P 

t - capt 


A P 


del 


(2-1) 


where 


p capt = P ower incident on and captured by antenna 
P del “ power delivered to load 

Since the received or captured power includes both incidence and 
re-radiation effects r it follows that 


P capt P del + P J 


( 2 - 2 ) 


where 


Thus 


Pj = power lost as Joule heating 


L a = 1 + p 


del 


so that 1 <_ L a < 


( 2 - 3 ) 

( 2 - 4 ) 


A radome inserted between an antenna and an incident field will 
introduce an insertion loss 


l r p 


B 


( 2 - 5 ) 
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where 


p = power received before radome insertion 
P B — power received after radome insertion 
It is observed that P > P_ if the radome does not interact with 

X3 A 

the near-field distribution of the antenna , since radome dielectric 
materials have non-zero a.c. conductivities. Thus, 


1 <_ ® 


( 2 - 6 ) 


Since antennas and radomes are lossy networks , it is clear 
that they will introduce a small amount of their own noise power, 
even when no external wave is incident. Generalizing to a two- 
port network, as in Fig. 2-1, we may solve for the apparent 
antenna brightness temperature, i.e. 


Emission Absorption 


T 1 
B 


1 A C 

» * x> t a + r 


(2-7) 


where it is assumed that an unpolarized random noise signal of 
equivalent temperature T g is incident on an isothermal antenna 
at thermometric temperature T A and of loss L. The first term in 
(2-7) is the component of the apparent brightness temperature 
due to the self-generated noise in the antenna and the second 
term is due to the absorption attenuated equivalent temperature, 
or received noise power. If the incident brightness temperature 
distribution is T sky (9,<f>) and the antenna gain pattern is G(e,<|>), 
then the second term of (2-7) is 



4 IT 


T sky (e,<f>) G (0 , $ ) dR 


(2-8) 
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2.2 Directivity/Gain Technique 

It has been unofficially estimated by NBS that both the 
gain and directivity of moderate to high gain antennas can be 
measured, using the near-field sampling technique [Kerns, 1970], 
to an accuracy of 0.1 dB (at the 3 0 point) throughout the 
frequency range L to K -Band. This is assumed to be the state of 
the art. We wish to examine the effect of such gain and directivity 
measurement errors on the radiometer antenna temperature 
uncertainty . 


The power gain of an antenna is related to its directivity 
by 



where 


(2-9) 


G = power gain (dimensionless) 

D = directivity (dimensionless) 

= antenna loss 

We next calculate the effect of systematic and independent 
errors AG, AD on the insertion loss error AL*. From (2-9) 


L, 


( 2 - 10 ) 


Differentiating and using the worst case formula. 


AL_ = —AD + ^-AG 
AG q 2 


It is assumed that there is zero systematic error. 


( 2 - 11 ) 
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which can be written as 


AL A = | < AD + LAG > 


( 2 - 12 ) 


It is next assumed that G and D are independent random variables, 

normally distributed with zero mean so that L - L a . Letting 
2 2 2 a - 

0 h r a G* a D be t ^ ie var i ances of L, G and D respectively, we 

can determine air by the familiar quadrature formula, i.e„, 

Jj 



’ 3L a’ 

2 

-2 

3L A 

SG 

ff G + 

. 30 



(2-1 3 ) 


which yields 


a 


L 


1 

G 





2 2 
L A a G 


(2-14) 


Since the errors are normally distributed we can express the 3 a 
error as 


3o L ■ I J a D + L A og2 


(2-15) 


The incremental errors o_ , a„, and a^ can be related to their 

ii b D 

corresponding decibel errors by 


3a. 


6L(dB) = 10 log (L + 3a L ) -10 log L = 10 logll + -j~ 


(2-16) 


If it is assumed that Sip(dB) ~ SG(dB), then it can be shown that 


6L(dB) = 10 log 


<5D 


i + 'xfi’uo 10 -i) 


(2-17) 


When SD - 6G = 0.1 dB, (9a) yields 


6L = 0.14 dB 


(2-18) 



11 


For example, when L & = 1 dB, 3 <j l = 0.041. 

We next consider the effect of the insertion loss error on 
the apparent antenna temperature, T R . The noise temperature T R 
at the terminals of an antenna is composed of two parts: (1) the 

integrated brightness temperature T g , diminished by the antenna 
loss and (2) the emissive temperature of the lossy antenna structure 
itself. 


Egn. (2-7) can be written as 


B 



T, 


(2-19) 


Assuming that T g and are known, the effect of a random error 
ct l on T g is given by 



< T A - 


T=> -T 

U * 

l a 


( 2 - 20 ) 


Assuming that the insertion loss uncertainty is entirely caused 
by the random SG = 6D = 0.1 errors (typifying the state of the 
art), eqns. (2-16), (2-17) and (2-20) can be used to calculate 
o T ^. For example, a 1 dB insertion loss (L^ - 1.259} corresponds 
to 5L = 0.14 dB (3a L = 0.041). Assuming T g = 10 K* and = 300 K, 


a- = (300 - 10.0) °*Q 41 7 = 7.5 K (2-21) 

i B 1.259 

For antennas with low insertion loss such as scalar horns, - 1 

A 

and a 0.14 dB insertion loss error would correspond to - 1 + 0.Q32 


This would be a typical figure for the antenna looking at the cold 
sky. 
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However, since >. 1, the effective 3 a L = 0< 5) 32 = 0,016, cor- 

responding to a_ = 4,6 K, Using this approach, we obtain the 

B 

curve of Fig. 2-2 which plots the antenna temperature uncertainty 
vs. the insertion loss for two assumed values of G, D measurement 
error. 


The preceding results demonstrate that in the calibration 
procedures where the measured insertion loss is used to deduce 
the integrated brightness temperature, relatively small errors 
in insertion loss can cause unacceptably high antenna temperature 
errors. In particular, when the insertion loss is inferred 
from the D/G ratio, independent random errors of 0.1 dB in the 
measurement of D and G correspond to temperature errors of 
5 - 9 K for typical values of insertion loss. 

Thus, single-frequency direct measurement techniques are 
not yet sufficiently accurate for the determination of the 
insertion loss in a microwave radiometric calibration situation. 
In addition, this method suffers from the disadvantage that an 
extremely large near-field sampling device would be required for 
large radome-covered radiometer systems such as PMIS or MFMR, 
and that to date no such large system has been built. 

2 . 3 Cryoload Technique 

A second technique which is usable for smaller antennas 
such as horns makes use of a cooled microwave absorbing 
hohlraum or cooled blackbody enclosure, as shown in Fig. 2-3. 

If the horn views an isothermal brightness distribution T_ in 

S 

a perfectly absorbing (and emitting) medium, then solving for 
from (2-7). 


INSERTION LOSS (dBJ 


Fig. 2-2. Apparent temperature error (AT„) versus antenna loss (L, ) 
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where 

= antenna loss 

T_ = thermometric temp, of antenna ^K) 
(assumed here to be isothermal) 

T g = brightness temp, of cryoload (K) 

Tg = uncorrected brightness temperature 
measured by radiometer (K) 


( 2 - 22 ) 


Blume and Swift [1972] of NASA LRC have used this technique to 
calibrate the loss of an S-Band horn and have achieved an 
absolute accuracy of + 1 K. 


The technique has the advantage of being inexpensive and 
accurate for relatively small antennas (diameters less than 6'), 
but as yet no isothermal cryoloads have been developed which are 
large enough to accommodate such large radiometer antennas/radome 

3 

systems as PMIS or MFMR where 30 m volumes are encountered. 


2.4 Bucket Technique 
2.4.1 Theory 

It is clear from (2-22) that the antenna thermometric 
temperature (T^) must differ substantially from the integrated 
source brightness temperature (T g ) in order for to be measured 
with acceptable accuracies. The bucket technique achieves this 
by placing the antenna in a large reflecting enclosure which 
blocks thermal emission from surrounding terrain, thus allowing 
the antenna to receive atmospheric noise only. The equivalent 
sky temperature can be calculated from radiosonde data by 
calculating the radiative transfer of an assumed 2.7 K cosmic 
background temperature through a clear atmosphere with both 
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water and oxygen resonant molecular constituents [Paris, 1971] . 
This technique is shown in Fig. 2-4,. 


In general. 


T = 
S 


JJI 

Af 477 


T sky (e,<j>,f) d ( 0 , <f> , f ) F(f) an df 


(2-23) 


4tt F(f) df 




where 


T sky 

F 

Af 


integrated antenna brighness temp. (K) 
sky brightness temp, distribution (K) 
spectral transfer function of radiometer 
system (Hz 
bandwidth (Hz) 


If the bandwidth is much smaller than the pressure-broadened 
spectral widths of the ^0 and (^ emission lines and if the 
antenna pattern is independent of frequency across the bandwidth, 
then (2-23) reduces to 

T S ” I, [[ T sky (9 ' ,,, ' f o ) D "•'♦'V ds! (2 - 24) 

4tt 

where f is the center frequency. If the sky temperature 
distribution is essentially constant across the main beam and 
first few sidelobes of the antenna in the bucket, then (2-24) 
reduces to 



' f o> 


(2-25) 


which is the form used in all subsequent calculations 
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2.4.2 Errors due to Bucket Emissivity 

It is apparent that (2-25) implies that the antenna in the 
bucket is assumed to be surrounded by an isothermal brightness 
temperature distribution, i.e., the bucket has a net zero emis- 
sivity. Also, it is tacitly assumed that the bucket is large 
enough so that there is no mutual coupling between the test 
antenna and the bucket. Assuming an emissivity of .001 for 
aluminum, the bucket brightness temperature would be [Carver ,1973] 

T bk = e T & = (.001) (290) = .29 K (17) 

Assuming that the main beam views a 10 K sky and a 95% bucket 
efficiency*, the contribution from the sky plus bucket would be 


T b n T sky + (1 “ ( T BK + T sky^ 


T sky + (1 " t BK 


= 10 + (.05) (.29) 


= 10.01 K 


Thus the error in neglecting the emissive temperature of the 
bucket walls is only .01 K. Even under severe oxidation conditions 
in which the emissivity might increase by a factor of 10, the 
error is still less than 0.2 K. 


The bucket efficiency n is the percent power received by the 
antenna which is not reflected by the bucket. 


2.4.3 Bucket. Shape Criteria 
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How large does the bucket need to be in order to avoid 
first order mutual coupling effects? The answer to this depends 
on the near-field distribution of the worst-case antenna (s) being 
used, i.e., the antenna (s) having the strongest fields near the 
bucket walls. For the case of PMIS/MFMR, the L-Band array and 
the X-Band array are electrically closest. These near-field 
distributions can be estimated using Hansen's [1964] calculated 
curves for uniform and tapered distribution horns, and setting 
a criterion that the distance from the antenna to any bucket 
wall should be large enough so that the near field reflected 
power is at least 30 dB below the main beam on-axis power at 
the same distance. 

2 v 4 4 Bucket Effect on Antenna Pattern - An X-Band Experiment 


A small model bucket was built using this approach, scaled 
for use with an X-Band standard gain horn, and with dimensions 
shown in F3g. 2-3. The power patterns of the horn in free space 
and inside the bucket are compared in Fig. 2-6 where it is seen 
that the effect of the bucket is to eliminate sidelobes beyond 
the bucket shadow boundaries and to introduce a distortion to 
the free-space pattern. The horn is relatively large in 
comparison to the bucket so that the perturbation is most likely 
due to the phase interference of multiply edge-diffracted rays 
and the unreflected incident rays. 

2.4,5 Description of Full Scale Bucket 

The full-scale bucket was made much larger than any antenna 
to be put in it, so that no appreciable perturbations in the 
lit- zone pattern is expected. Dimensions are given in Fig. 2-7 
along with sketches of the antenna positioner and PM1S/MFMR 
systems. 
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Fig. 2- 



5. Scale model bucket for X-Band standard gain horn. 
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Fig* 2-6. Standard gain horn antenna pattern. 


j 





23 


The bucket is an aluminum foil covered wooden truncated and 
inverted pyramid, with a reinforced concrete base pad. It is 
located about 9 km east of the NMSU campus at an altitude of 
1.46 km on "A” Mountain. The photographs of Fig. 2-8 - 2-11 show 
phases of the bucket construction. Equipment access is through 
a 3 m x 6 m door, as shown in Fig. 2-12, with special lifting 
rigs available for the PMIS and MFMR antenna and radome systems, 
as shown in Figs, 2-13 - 2-16. Personnel access is through a 
small removable door on the south side. Antenna/radome mockups 
are mounted on a Scientific-Atlanta as over el positioner affixed 
to the top of a steel support tower, with center-of -gravity 
maintained on the tower axis by use of counterbalance assemblies. 

Electronic and control equipment is housed in a large 
temperature controlled concrete block building (see Fig. 2-17) 
about 15 m south of the bucket with an interconnecting cableway. 
During operations, a weather station is used to give temperature, 
relative humidity, pressure, wind speed and wind direction. An 
all-weather gravel road permits truck transport of equipment to 
the facility. 



2-8 Pre-fabrication and Painting Bucket Parts 
(12 Nov. 1974) 


2-9 Cementing Aluminum Foil to 
Bucket Interior Surface 

(13 Nov. 1974) 


ORIGINAL- PAGE $ 
OF POOR QUALITY 







Fig. 2-10 Bucket Assembly Underway 
(14 Nov. 1975) 


2-11 South Face of Bucket Showing Braces 
and Rigging to Sustain 100 mph Wind 

(Jan. 1975) 
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Fig. 2-12 PMIS in Position for Test - 
Main Bucket Door Open 

(Jan. 1975) 



Fig. 2-13 PMIS in Mockup on Sling and Handling Dolly 

(6 Feb. 1975) 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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Fig. 2-14 MFMR Radome in Handling Sling 
(20 Feb. 1975) 
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Fig. 2-15 




Fig. 2-16 Mounting PMIS - Note Counter Weight 

(24 Jan. 1975} 
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Fig. 2-17 PMIS Equipment Setup 
(Dec. 1974) 



Fig. 2 18 flFMR with Absorber Mounted for Test 


original pacp to 
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(Feb. 1975) 
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CHAPTER III 


ANTENNA/RADOME DATA REDUCTION TECHNIQUES 
3 . 0 INTRODUCTION 


This chapter discusses the theoretical framework for the 
inference of antenna and radome loss as well as the practical 
techniques used for the PMIS and MFMR systems. 

It is worth repeating that the loss values and L R are 
used in establishing the absolute accuracy of an integrated 
radiometer (antenna plus receiver) and that the need for these 
values must be established by the user. How accurate must the 
inferred brightness temperatures be? At present, accuracies of 
+ 1 K are considered to be very good with receiver sensitivities 
between 0.1 K and 1.0 K achievable using integration times rang- 
ing from 100 - 1000 milliseconds. But can a user establish a 
need for a + 1 K accuracy? Will + 3 K suffice? This is a dif- 
ficult question, having to do with the adequacy of radiative 
transfer models, the availability of sufficient ground truth 
in controlled areas, aircraft/spacecraft platform stability, etc. 

It is clear, however, that relative loss values are extremely 
important. In the case of PMIS, there are 44 beam positions, 
electronically scanned, and it is necessary to establish the 
loss value of a given beam position relative to its adjacent beam 
positions with maximum precision. For MFMR, the variation in 
loss with pitch angle change is of greater immediate importance 
than the absolute value of those losses. 

3. 1 Theoretical Foundation for Loss Calculations 

A calibrated radiometer viewing an isothermal source T g 

indicates an uncorrected brightness temperature T * which is 

B 
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greater than T g , due to the added emissive contributions of the 
intervening waveguide, antenna, and radome. Considering first 
the radiative transfer when there is no radome (Pig. 3-la) , we 
have [Kraus, 1966] 


T t 

B 


(1 



+ 


1_ 


emissive 
contr. from 
waveguide 



antenna 

absorption 

term 



k ■ j — ■ J 


emissive 
contr. from 
antenna 


(3-1) 


where 


T' = uncorrected brightness temp, measured 

D 

by radiometer (K) 

Tg - sky brightness temp. (K) 

T a = antenna kinetic temp. (K) 

T w - waveguide kinetic temp. (K) 

L a = antenna loss (L^ > 1) 

= waveguide loss (L^ > 1) 


It is assumed that the waveguide loss is known from 
laboratory measurements. Solving for L^, 


r 


t a “ T s 

T* + (L r . - 1 ) T,. - L r ,T ’ 


w 




W B 


(3-2) 


Turning now to the radome loss L R (Fig. 3-lb) , the radiative 
transfer equation is 


1 1 

T* = (1 - =—) T + =— 

B 1 V » h, 


a - i-)T + f- 

L a A l a 


1 s 

(1 - =-)T „ - 

L R R L R J 


(3-3) 







Solving for L 


t r + t a' l a - 11 + Wh, - » - We 


The coupled equations (3-2) and (3-4) form the basis for the 
measurement of the antenna loss and radome loss. It is assumed 
that the L A value in (3-4) is the same as that calculated in 
(3-2), i.e., there is no mutual coupling between radome and 
antenna • 

The thermometric temperatures of the antenna, radome and 
waveguide structures are read by imbedded thermistors. The sky 
brightness temperature is calculated by techniques described in 
Chap. IV and the uncorrected brightness temperature T* is that 

iJ 

indicated by the calibrated radiometer. 

In practice, both random and systematic errors are introduced 
into (3-2) and (3-4) with concomitant effects on both precision 
and accuracy. These are discussed in detail in Chaps. V and VI. 

3 . 2 PMIS/MFMR Data Systems 


Both PMIS and MFMR are operated on a Lockheed P-3A aircraft, 
along with other remote sensors such as multi-spectral optical 
scanners, etc., with all electronic data being recorded on 
magnetic tape in PCM format. Thus, the output of both PMIS and 
MFMR is in PCM counts with integration times of 100 ms being 
typical during flight operations. However, an integration time 
of 1 minute was used during the loss calibration tests so that 
an improved sensitivity would be obtained. 




The sensitivity of the radiometer is continuously monitored 
by the use of two internally generated equivalent temperatures 
known as the calibrate and baseline modes, with corresponding 
PCM counts and C , respectively. Assuming the radiometer is 

V- D 

linear, the uncorrected brightness temperature is given by 

(3-5) 

where 

T^ & AT = constants to be determined (K) 

C A = one minute average of data counts 

C_. - one minute average of baseline counts 

C c = one minute average of calibrate counts 

as shown in Fig. 3-2. 

T^ and AT are determined by connecting known noise temperatures 
to the radiometer input port and by computing the best linear 
regression fit to the resulting data, or if only a hot and cold 
load are available, by solving two equations and two unknowns. 

Figs. 3-3 and 3-4 show the calibration history of T^ and AT for 
the PMXS vertically and horizontally polarized channels, over 
the Jan. 20 - Feb. 6, 1975 time period. The horizontal channel 
shows much more long-term instability than the vertical, which 
was the basis for recurring difficulties in data reduction. 

Fig. 3-5 shows a similar history of T^ and AT for MFMR, comparing 
the initial calibration on Feb. 8, 1975 to the final calibration 
on Feb. 21, 1975. The K-Band (22.05 GHz) receiver was not working 
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on. Feb. 21, so that only an initial set of values are shown. The 
apparent downward drift in the T^ value for K u -Band was caused 
by not allowing sufficient warmup time before calibration. 

In addition to the PCM counts corresponding to the radio- 
metric data, both PMIS and MFMR software require the monitoring 
of other housekeeping data, such as antenna, radome and waveguide 
thermistors, AGC voltages, electronics enclosure temperatures, 
etc. These data are multiplexed and recorded on magnetic tape 
for subsequent processing, as shown in Fig. 3-6. Data tapes 
recorded at the bucket site are brought to the computer facility 
at the PSL building on campus. Data on the tapes are then decom- 
mutated, and transferred to storage addresses within the IBM 
370 computer, all of this being under control of the IBM 
SYSTEM 7 process controler. 

Data reduction then proceeds under three phases of software: 
PMIS/MFMR -1, 2, and 3 as shown in Fig. 3-6. Phase 1 strips and 
sorts data from the tape, calculates means and standard deviations 
and prints significant data in counts, as shown in Figs. 3-7, 3-8. 
Phase 2 converts radiometric housekeeping data into appropriate 
engineering units, as shown in Figs. 3-9, 3-10. Phase III then 
calculates loss values and associated standard deviations by 
using eqns. (3-2) and (3-4), and prints the results as shown in 
Figs, 3-11, 3-12. 

The detailed instructions for the reduction of PMIS and MFMR 
data tapes are given in the Appendix. The software written for 
this purpose is extremely complex and powerful. However, it must 
be complemented by a parallel intervention of the project engineer 
and programmer in order to monitor the health of the various 
components of the PMIS/MFMR systems and to properly intermesh the 
values (eqn. 3-2) into the computation of L R (eqn. 3-4) . 
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Fig. 3-7 PMIS-1 Sample Output . 
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Fig. 3-8 MFMR-1 Sample Output. 


tfs- 

u> 



44 


P 1SK RAW DESCP IPT ION CHANNEL DISK ENGINEER | NG UNITS FORMAT IDP OR EOP 


ADDRFSS 
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NUMBER 
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SUBSCR 

508000 
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BEAN POSITION 

3- 1-0 

508050 

22 



I 

1 

508001 

66 

V/H PATIO 

3- 2-0 
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I 

2 

508002 
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F 

3 
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0 
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F 

4 
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I 

5 
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STO vertical 

3-22-D 
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6 
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0 
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F 
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F 

12 
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c 

F 

13 
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F 
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F 
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F 
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F 
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c 

F 

20 
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0 
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c 

F 

21 
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246989 
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c 

F 

22 
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c 

F 

23 

508023 
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c 

F 

24 

508024 
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HORIZ ENCLOSURE TEMP 

3-58-A 
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OEG 

c 

F 

25 
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3-89— A 
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F 

26 
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c 

F 

27 
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Fig. 3-9 PMIS-2 Sample Output 



45 


LOGIP = MOtn VALID DATA cone = 000007.45 

START TIM8 = 2 HPU8S 45 WTNIJTES 46.03 SFCnNDS END M MC = 2 HOURS 46 MINUTES 46.16 SFCCNDS 

RAW TOTAL DATA FOB THIS FILE STARTS AT DISK ADDRESS 247S0C’ 

r Nr.lN=FKTNC PA T A FUR THIS F»Lf SJ ARTS AT DISK ADDRESS 2475EC 


• 

°AW 

. SFNSCP SYSTEM 

. THAN 

. 'NGINFFn fdjo 

. 
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10 

• 
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• 

R. /5P0F 

02 
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li 

* 
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KU SWITCH TFMp 
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* 
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C2 

DFG K 

12 

* 
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0? 

DEG K 

13 

• 
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* 
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02 

DEG K 

14 

• 

135028 
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2A-9-A 

* 

2 . 9969 r 

00 

VOLTS 

15 

• 

12955A 


KA-BAND COLO REF. TFMP. 


2 A— 1 5- A 

* 
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02 

DEG K 

16 

• 

13A860 


KA HOT P f F . TEMP. 


2A-16-A 

• 

7.R207F 

P3 

DEG K 

17 

« 

121063 


K1 ANTFNNA TFMP. 


2 A-17-A 

• 

7.7736F 

02 

DFG K 

15 

* 

1295P0 


KA SWITCH TEMP. 


2 A-lfi-A 

« 

3.735TF 

02 

DFG K 

19 

• 

118685 


KA CH. 1 WAVE GUIOF TEMP. 


2 A- 19- A 

* 

2.7774= 

0? 

DEO K 

20 

• 

12008® 


KA CH.2 WAVE GUIC TFMP. 
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* 
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• 
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• 
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* 
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K SWITCH TEMP, 
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■ 

7.7579F 

07 

DEG K 

27 

* 

1 1 9034 


K CH.2 WAVE GUIDF TEMP 


2A-22-A 

■ 

2.7579F 

0? 

DFG K 

28 

• 

127653 


K-BANp AGC VOlTf.GF 


2 A- 33- A 

* 

2.0; ?? f 

00 

vm. T? 

29 

• 

135374 


L-3AND CDLP RFF. TFMP. 


2A-39-A 

■ 

2.70S3F 

07 

DEG K 

10 

* 

125124 


L HOT RFF. TFMP, 


? A— A 0- A 


3.B129F 

07 

nrr. k 

31 

• 

116456 


L ANTFNNA TFMP. 


2 A— 4 1 — A 

» 

2.7810F 

02 

nee, k 

32 

• 

127600 


L SWITCH TF«p. 


2 A- A?- A 

• 

1.21B9B 

07 

DEG K 

33 

* 

107214 


L WAVE GUIDF TFMP. 


2 A-A’-A 

* 

2.7973F 

0? 

DFG K 

34 

B 

116775 


L-RAND AGC VHLTAGF 


2A-A5-A 

• 

2.5909E 

00 

volts 

35 

m 

130243 


L-PAND FILTFP TFMP. 


2A-A6-A 

C 

1.2423E 

C? 

DEG K 

36 

* 

n 


RADOME THFPWIST'”? N'T.l 


29- 51 - A 

* 

2.1800F 

02 

OTG K 

37 

• 

o 


NO, A 


7B-54-A 

• 

7. IROOF 

02 

n FG K 

3H 

a 

0 


NO. 5 


28- 65- A 

* 

2 . i floor 

07 

DEG K 

39 

4 

0 


MI. 6 


2 a — 54- A 

■ 

7.18C0F 

07 

ntG K 

AO 

a 

7 


PAOOMF THERMISTOR NO. 7 


2B-57-A 

* 

7.1H00F 

02 

DEG K 

Al 

• 

5 72 00 


Mnnc r pilap i?at ton switch 


2B-l-n 

■ 

7.60001 

O'* 

rrjnr 

A2 

* 

2 5844 


ANTENNA FLEVATION ANGl F 


28-er-A 

* 

1.59728 

07 

AVG DFG 

A3 

• 

768652 


SIGMA Ktl— BAND CH.l 


2 A— 1 - A 

• 

1.25178 

01 

COUNTS 

AA 

• 

763005 


SIGMA KU-flANn CH.2 


2A-2-A 

* 

2.61 IAE 

00 

COUNTS 

A5 

* 

628911 


SIGMft KA-BAND CH.l 


ZS-17-A 

* 

1 .f S57 c 

PD 

COUNTS 

A 6 

• 

623734 


SJOMA KA-8AN0 CH.2 


2 A- 1 A-A 

* 

2 . 5 1 45 c 

no 

COUNTS 

A 7 

V 

734540 


SIGMA K-BANO CH.l 


2A-25-A 

* 

6.34028 

DO 

COUNTS 

AH 

* 

826405 


SIGMA K-0ANO CH.2 


2A-26-A 

■ 

2.5567F 

CO 

COUNTS 

A 9 

• 

605619 


SIGMA L RAND 


2 A— 3 7- A 

• 

1.62831 

01 

COUNTS 

30 


220 









51 

• 

581 









5? 

* 

«003 









5? 

« 

170061500 









5 A 

• 

1711 862 n 0 










Fig. 3-10 MFMR-2 Sample Output . 
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Fig. 3-11 PMIS-3 Sample Output. 


CD 


ij 




1 nr. m 

= *013 

VALT9 

9ATA C90F 

= 00000255 



STAOT 

= 

? HOU-5 

55 HlMJTfS 

56.02 SECONDS 

c NO time = 2 HQU9S 56 “tMUTCS 16.06 SFCCAIPS 



ANTENNA LOSS 





$KV 

T9 

SIGMA 

ANT. 

SIGMA 

ANT. 


TCMp 


T R 

LOSS 

LA 

LOSS 


tK J 

<KJ 

{ Kt 

LA 


10*1 

K V'- 1 

6. 17 

57.65 

5.9100 

1.1132 

0.0272 

0.4655 

K'J-2 

6.17 

39.05 

5.8059 

1.1C69 

0.027? 

0.1512 

KA- l 

10. °1 

59.93 

5.9572 

1.1210 

3.0256 

0.5962 

K A- 2 

10.11 

5B.47 

5.1687 

1.1155 

0.0263 

C .4703 

K-t 

12.16 

45.48 

5.7651 

1.0826 

0.0269 

o.i4‘»a 

K-2 

13.16 

53.27 

5.5761 

1.0777 

0.0257 

0.3251 

I 

4.?9 

90.56 

5.2719 

1.2661 

0.0180 

1.3548 


Fig. 3-X2 MFMR-3 Sample Output . 


48 


CHAPTER IV 


ESTIMATION OF SKY BRIGHTNESS TEMPERATURES 
4 . 0 INTRODUCTION 


Eqn. (2-23) and its special case (2-25) form the basis for 
estimating the isothermal brightness temperature seen by the 
antenna. The purpose of the bucket is to equate this temperature 
to the sky brightness temperature at zenith. 

The calculation of this temperature using the spherical shell 
model of Paris [1971] depends on the availability of contemporary 
radiosonde data at the bucket site in clear sky conditions, and 
it assumes negligible random fluctuations (within-the-hour ) . In 
practice, however, logistic considerations make it difficult to 
obtain soundings at the exact desired time and/or place. These 
and other factors can introduce systematic errors in which 

become greatest at the water resonance frequency of 22.235 GHz. 

This chapter discusses the general technique and pitfalls 
involved in this approach. 

4 . 1 Theoretical Background 


The low-level noise power emanating from the atmosphere in 
the microwave spectrum results primarily from absorption and 
re-radiation by water and oxygen molecular constituents with an 
extremely weak cosmic background ( T cosni j_ c = 2 *7 K) , believed by 
many cosmologists to be due to radiation from the remnants of a 
diffuse expanding primiordal fireball. Superimposed on this are 
occasional localized radio objects such as the sun, the galactic 
center and various point sources. Almost all of these, however, 
have flux densities which decrease rapidly with increasing 
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frequency and may be safely ignored at X-Band and above. At 
1.4 GHz it is possible to detect solar radiation with only mod- 
erate gain antennas, so that mid-day calibration measurements 
should be avoided at L-Band. 

For PMIS and MFMR, the task of primary concern is the 
modelling of radiative transfer within the atmosphere. The 
model constructed by Paris [1971] uses concentric spherical 
shells to describe the radiative transfer, as shown in Fig. 4-1. 
The equation of radiative transfer through one shell (shown in 
inset) is 


T 


BAZ 


m -otAZ sec 9 . -a AZ sec 0. 

= T BQ e + T(1 - e ) 


where 


(4-1) 


T 

T 

T 

a 


Bo 


BAZ 


AZ 

6 


incident brightness temperature (K) 
mean thermodynamic temperature (K) 
transmitted brightness temperature (K) 
volume absorption coefficient (nT 1 ) 
thickness of shell (m) 
angle from zenith (rad) 


By summing over a sufficient number of shells (with thickness 
equal to or less than 30 mb pressure) , the total incident bright- 
ness temperature can be computed for any station altitude or 
pressure. 


The volume absorption coefficient a is due to absorption 
by both water and oxygen and is modified to include pressure 
broadening effects using the model of Van Vleck and Weisskopf 
(1945) . The necessary input data for (4-1) are taken from radio 
sonde profiles of air temperature (T) , pressure (p) and relative 
humidity (R.H) , which are then converted to more useful forms. 



cosmic 


Typical shell 



It is necessary to include sounding data to 10 mb in order to 
accurately assess the sky temperature at 22 GHz where total 
precipitable water is of great importance. 
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Unfortunately, however, many soundings are terminated at 
300 mb so that it is necessary to simulate profiles above this 
level. This simulated tropopause effect is shown in Table 4-1. 



It has been found through computer simulation that these numbers 
are not critical, even at 22 GHz, but that some reasonable simu- 
lation must be made. 

4.2 Local Topography 

The radiometer calibration facility is located approximately 
7 km east of the NMSU campus in Las Cruces, and is at an altitude 
of 4816* (1.468 km), as shown .in Fig. 4-2. Most of the regular 
soundings (daily, at 0200 MST) are taken at White Sands Missile 
Range WSD site, although other soundings are taken occasionally 
at the SMR, LC-36 and Airport sites shown. Between the WSD site 
and the A-Mtn. bucket site is the Organ Mountain range with some 
peaks rising to 9000' (2.743 km). 


NEW 

MEXICO 



Fig. 4-2. Location of NMSU/PSL 
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Can the regular soundings at WSD be used to adequately 
estimate the sky temperature at the bucket site? The answer is 
a qualified yes, with the proviso that clear sky conditions be 
prevailing, that the WSD sounding be within 3 hours of the 
measurement , and that the wind be predominantly westerly. This 
technique has evolved as the result of systematic study of the 
correlation between sounding data in the area, and of the cor- 
relation between total precipi table water (W ) and the sky 

r 

temperature T g ^y. For example. Fig. 4-3 compares contemporary 
soundings at A-Mtn. and WSMR at 0200 MST on Feb. 12, 1975, and 
shows a higher temperature and slightly higher humidity on the 
WSMR side of the mountains. The effect of this difference on 
the T slc y spectrum is shown in Fig. 4-4, where the sensitivity to 
total precipitable water at K-Band becomes immediately apparent. 


These data are referenced to an average A-Mtn. surface pressure 
of 850 mb. The effect of smaller pressures (higher altitudes) 
is shown in Fig. 4-5 in the T s ^ spectrum. As the pressure 
broadening decreases, water and oxygen line shapes become more 
distinct, and the 2.7 K asymptotic value is approached at lower 
frequencies. 


The correlation between the zenith sky temperature and 

total precipitable water (W ) is shown in Fig. 4-6 for 18, 22 

P 

and 37 GHz. The total precipitable water is found by integrating 
from the atmospheric top (10 mb) to the station pressure of 
850 mb, using sounding data and standard meteorological expressions. 
At 18 and 37 GHz the computed zenith sky temperature is closely 
correlated to the indicated linear fits with slopes of 0.28° and 
0.38° per mm respectively. At 22 GHz the slope of 2®/mm is much 
steeper and the correlation is not as good, with one point depart- 
ing 2.5 K from the standard curve. The points shown on this scat- 
tergram were computed from sounding data in the winter time period 
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Fig. 4-3. 


Radiosonde profiles - Feb 
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Pig. 4-5. Zenith sky temperature versus frequency. 
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Nov. 1, 1974 to Feb. 20, 1975. A histogram showing the distri- 
bution of 22.05 GHz temperatures is shown in Fig. 4-7 and 
indicates the most recurrent temperatures to be in the 14-17 K 
range [Carver, Cooper and Paris, 1975] . 

Even under clear sky conditions, the zenith sky temperature 
can vary widely over a period of several days, as shown in Fig. 
4-8, which compares A-Mtn. temperatures to WSMR (WSD) values in 
the Feb. 12-14, 1975 time period. Excursions over 10 K are 
evident at 22 GHz (WSD) in a little less than two hours. 

Since 22 GHz is the most difficult frequency from the stand- 
point of prediction, a key question is whether the zenith sky 
temperature can be estimated with sufficient accuracy using only 
total precipitable water (W^) as a basis. Fogarty [1975] has 
very recently reported on the correlation between 22.2 GHz 
measured zenith temperatures and surface dew point temperatures 
over the Oct. 1972 - Nov. 1974 time period for a Brazilian 
coastal zone site 850m above MSL. He observed large seasonal, 
daily and even hourly variations in the zenith attenuation 
with typical clear sky values of 0.64 dB in the winter and 
1.70 dB in the summer, corresponding to 40 K - 88 K temperature 
variations respectively. He concludes that, contrary to previous 
recommendations [Sullivan, 1971] , the prediction of zenith sky 
temperature by surface dew point temperature or total precipitable 
water is much less accurate than by use of the "tipping" method, 
wherein the sky temperature is observed at various angles and 
calculated by slope techniques. The difficulty with Fogarty's 
assertion is that no explanations are offered to explain how he 
lculates total precipitable water, although it is apparent that 
c is not the total integrated used in the Paris approach. 


ZENITH SKY 


Fig. 4-7. Histogram of 










ZENfTH SKY TEMPERATURE (K) 
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It is believed that the proper use of radiosonde profiles 

with an adequate radiative transfer model leads to an acceptable 

correlation between the zenith sky brightness temperature and W , 

I? 

as indicated in Fig. 4-6. Clearly, a direct calculation of 
using the SKYTEMP program written by Paris with timely sounding 
data from A-Mtn offers the most accurate means of estimation. 
However, on-site soundings are relatively expensive and logisti- 
cally difficult so that it is desirable to use the regular 0200 
WSD soundings where possible. 


By comparing the integrated W values from several coordinated 

]? 

soundings from Las Cruces Airport, Small Missile Range, and 
Launch Complex - 36, it has been found that for clear sky condi- 
tions and prevailing westerly winds, the value is about 1/2 mm 
higher on the east side of the Organ Mountains {see Fig. 4-2) , 
so that from Fig. 4-6, T sk y should be higher on the west side of 
the mountains. This is confirmed in Fig. 4-8 where A-Mtn. and 
WSD temperatures are compared, with the brightness temperatures 
on the A-Mtn. side being consistently cooler. 


4 . 3 Variation with Zenith Angle 

Fig. 4-9 shows the variation in the computed T g j,y values at 
18, 37, and 22.05 GHz with angle from zenith. Both A-Mtn. and 
WSD curves (Feb. 12, 1975; 0200 MST) are shown. The greatest 
slopes are found at 22 GHz, although the variation is negligible 
over the 5° beamwidth of the scalar horn used. Thus the assumption 
that T g - ^g^ytO 0 )* as in eqn. 2-25, is a good one when a reflect- 
ing bucket enclosure is used to block emissive radiation from 
surrounding terrain. 


ANGLE 


Fig. 4-9. Comparison of computed 
from zenith. 



2/1 2/75, 0200 MST, A MTN 


2/1 2/75, 0200 MST, WSMR 








4.4 Errors in T. 


An extensive analysis of systematic errors in T^y using 
radiosonde data has been prepared by Paris (1975] especially 
for the PMIS/MFMR measurement program. His analysis assumes that 
the sounding is contemporary with the radiometric measurement 
and is taken at the A-Mtn. site and that random errors can be 
ignored. These error values are repeated in Table 4-2. 


Table 4-2 


Systematic Errors in T e 
(after Paris, 1975) 


Frequency (GHs) 


^Skv (K) 


1.4135 

10.69 

18.00 

22.05 

37.00 


0.27 

0.40 

0.82 

3.04 

1.84 
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CHAPTER V 


PMIS LOSS MEASUREMENTS 


5.0 INTRODUCTION 

The PMIS radiometer operates at 10.69 GHz as a dual- 
polarized imaging sensor using a planar phased array of crossed 
slots. The beam is scanned through 44 discrete positions on a 
conical surface and has a typical beamwidth of 2° corresponding 
to a gain of 35 dBi. Horizontally and vertically polarized 
components of incident radiation are processed by separate 
radiometers with beam switching and data initial data processing 
under control of a dedicated computer. In flight, integration 
times are variable and controlled by a feedback network. A 
radome covers the array when it is used on the P-3A aircraft. 

In the calibration testing phase, the radiometer’s intrinsic 
integration time was fixed (~ 120 ms) with an effective integration 
time of 1 minute being provided by computer processing of data 
tapes. The objective of the test was to provide antenna loss 
(L n ) values for all 44 beam positions for both vertical and 
horizontal channels, and radome loss (L R ) values for both 
radomes supplied, at each of the 44 positions. This totals 
to 264 separate loss measurements. Tests were run at night, 
normally from midnight to dawn, so that the regular 0200 MST 
radiosonde sounding from W3D could be used to calculate the sky 
brightness temperature. 

5.1 Mechanical Positioning Technique 


The PMIS array and associated radiometers were mounted in 
a plywood bomb-bay mockup which was in turn affixed to an 
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azimuth-over-elevation antenna positioner, as shown in Fig. 2-16. 
The positioner was set {according to a table in the operator 
logbook) so that electronic beam steering was compensated, giving 
a beam always pointed toward the zenith. 

A counterbalance assembly was used so that the center of 
gravity was nearly on the elevation axis. Angular setting errors 
of the positioner are negligible. 

5.2 Radiometer Receiver Calibration 


Both horizontal and vertical radiometer receivers were 
laboratory calibrated by the use of an X-Band waveguide variable 
temperature cold load, as shown in Fig. 5-1. A matched wave- 
guide load is immersed in a cryobath whose temperature can be 
varied from approximately 50 K - 250 K. This source of thermal 
noise causes an equivalent noise temperature T f to appear at the 
A-A* waveguide flange connected to the radiometer receiver. 
Thermistors {one is shown) are used to monitor temperatures 
along the waveguide. Corrections are made for waveguide 
absorption and internal emission, as well as mismatch corrections, 
in arriving at the A-A' flange temperature. 

This reference load has not to date been calibrated by 
NBS and it would be desirable to have this done. However, it is 
instructive to consider possible sources of error in order to 
arrive at a total systematic uncertainty in the flange temperature. 
These errors result primarily from inaccuracies in the cryobath 
temperature (+ .2 K ) , and flange mismatch (misalignment) errors 
(+ .5 K) with the resulting estimate of a + 1 K overall systematic 


error. 
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. Symbolic diagram of X-Band reference cold load for 
PMIS receiver calibration. 


Fig. 5-1 




In practice, the apparent flange temperature (equated 
to T') in 25 K steps from 100 K to 200 K, which is the normal 
operating range of the receiver. As the cryoflask mixing is 
changed, short-term local temperature gradients develop around 
the matched load so that T f requires several minutes (typically 
10 min.) to stabilize after making a 25 K change. An equation of 
the form 


'£i = T Bi " T 1 + l4T)X i 


(5-1) 


where 


X. - 
i 


C A - C B 


(5-2) 


C C - C B 


is used to solve for the calibration constants T^ and AT by 
using standard regression techniques to find the best straight 
line fit to (5-1) . T^ assumes approximate values 100 K, 125 K, 
150 K, 175 K and 200 K, and is read from the (827) address of 
the Varian computer. 


The T^ and AT "constants" varied slightly, due to receiver 
post-detection instabilities, as shown in Figs. 3-3 and 3-4. 

The horizontal receiver showed about 3 times the drift of the 
vertical receiver and generally provided a continuum of dif- 
ficulties. A further source of difficulty in receiver cali- 
bration was a high level of PCM noise, which was corrected only 
near the end of the measurement program. 


On Jan. 27, the horizontal receiver failed completely, so 
that all subsequent data from both horizontally and vertically 
polarized channels were taken using the vertical receiver front- 
end. 
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5.3 PMIS Data Flow 


Each of the 264 loss values were measured independently 
several times so that a reduction in the uncertainty was made 
possible by averaging. The number of independent measurement 
sets is shown in Table 5-1. 


Table 5-1 


Number of Independent Loss Measurement Sets 


LOSS 

Polarization 

Number of Sets (N g ) 

Antenna 

Vert. 

7 


Horiz . 

7 

Radome #1 

Vert. 

3 


Boris. 

4 

Radome #2 

Vert. 

3 


Horiz. 

4 

i 

1 — . 1 . • . 1 


The mean loss value for each set was then computed ac- 
cording to 



1_ 

N s 


N s 

l 

i=l 


M 


(5-3) 


where the superscript refers to the beam position number and 
N c is the number of sets. The standard deviation of the set was 
computed according to 


M 

0 


N S 


N „ 

l (l. m - l m ) 2 
i=l 1 


> 1/2 


1 


(5-4) 
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M 

It was found that the noise o was still too high using this 

procedure, particularly for the horizontal channel. To reduce 

the noise further, a three-point convolution process was used, 

wherein a loss value for beam M is averaged in a weighted sense 

with values of its two neighbors M-l, M+l. The weighting 

function is chosen as the standard deviation, so that L values 
• M 

having high (noisy) a values count less than those with more 

M 

repeatable values (low a ), i.e., the final quoted loss value 
is 


l m = 


M+l , -5J 

l (o N ) _ 1 L N 
M-l 


M+l 


I <a w > 

M-l 


N, -1 


(5-5) 


As an example, the antenna loss (horizontally polarized case) 
for beam position 4 is derived from measured data listed in 



t 
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Then 


L 


4 


1.557 , 1.554 
.04 .04 


704 + 704 + 


1.572 

.05 

1 

.05 


1.560 


This smoothing process essentially borrows statistical information 

from neighboring beams and is equivalent to convolution of 

M 

average loss values L with a three-beam weighted pulse function. 
5.4 PMIS Antenna and Radome Loss Values 


The final loss values, as derived by this technique, are 
tabulated in Table 5-3 and illustrated graphically in Figs. 5-2, 
5-3 and 5-4. 
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Table 5- 

‘A 






PMIS Loss 

Summary 




Beam 

Position 

Antenna 

Loss 

Radome 

1 Loss 

Radome 

2 Loss 


Vert. 

Horiz. 

Vert. 

Horiz . 

Vert. 

Horiz . 

1 

1.672 

1.555 

1.083 

1.130 

1.064 

1.125 

2 

1.677 

1.557 

1.083 

1.134 

1.063 

1.127 

3 

1.679 

1,558 

1.083 

1.136 

1.060 

1.129 

4 

1.679 

1.560 

1.082 

1.137 

1.059 

1.131 

5 

1.682 

1.572 

1.080 

1.138 

1.057 

1.133 

6 

1.685 

1.585 

1.079 

1.138 

1.056 

1.134 

7 

1.681 

1.585 

1.077 

1.139 

1.055 

1.135 

8 

1.677 

1.586 

1.077 

1.139 

1.055 

1.135 

9 

1.669 

1.589 

1.076 

1.139 

1.055 

1.134 

10 

1.667 

1.595 

1.074 

1.139 

1.054 

1.134 

11 

1.664 

1.604 

1.072 

1.139 

1.053 

1.132 

12 

1.666 

1.614 

1.069 

1.139 

1.052 

1.130 

13 

1.664 

1.634 

1.067 

1.138 

1.051 

1.128 

14 

1.661 

1.641 

1.066 

1.138 

1.050 

1.127 

15 

1.657 

1.642 

1.066 

1.138 

1.050 

1.125 

16 

1.655 

1.638 

1.066 

1.137 

1.050 

1.125 

17 

1.650 

1.642 

1.066 

1.137 

1.049 

1,124 

18 

1.646 

1.649 

1.066 

1.136 

1.048 

1.124 

19 

1.646 

1.646 

1.067 

1.134 

1.048 

1.123 

20 

1.645 

1.640 

1.067 

1.132 

1.048 

1.122 

21 

1.644 

1.638 

1.068 

1.130 

1.048 

1.121 

22 

1.641 

1.637 

1.070 

1.126 

1.048 

1.120 

23 

1.639 

1.645 

1.071 

1.123 

1.048 

1.118 

24 

1.639 

1.641 

1.071 

1.119 

1.049 

1.116 

25 

1.643 

1.656 

1.072 

1.116 

1.048 

1.114 

26 

1.644 

1.655 

1.071 

1.112 

1.048 

1.111 

27 

1.646 

1.659 

1.071 

1.109 

1.047 

1.108 

28 

1.649 

1.661 

1.069 

1.106 

1.047 

1.106 

29 

1.658 

1.660 

1.068 

1.104 

1.046 

1.105 

30 

1.662 

1.662 

1.066 

1.105 

1.046 

1.106 

31 

1.664 

1.651 

1.065 

1.110 

1.046 

1.108 

32 

1.665 

1.635 

1.065 

1.119 

1.046 

1.111 

33 

1.667 

1.624 

1.066 

1.125 

1.047 

1.113 

34 

1.672 

1.615 

1.067 

1.129 

1.048 

1.115 

35 

1.675 

1.604 

1.068 

1.131 

1.049 

1.117 

36 

1.679 

1.594 

1.071 

1.132 

1.050 

1.118 

37 

1.684 

1.592 

1.073 

1.134 

1.052 

1.119 

38 

1.680 

1.585 

1.075 

1.135 

1.054 

1.120 

39 

1.684 

1.579 

1.078 

1.136 

1.056 

1.121 

40 

1.688 

1.570 

1.080 

1.136 

1.057 

1.121 

41 

1.686 

1.562 

1.080 

1.135 

1.058 

1.122 

42 

1.681 

1.552 

1.079 

1.133 

1.059 

1.122 

43 

1.681 

1.552 

1.076 

1.126 

1.059 

1.122 

44 

1.6C1 

1.547 

1.075 

1.116 

1.059 

1.122 
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5. 5 FMIS Error Budget 
5.5.1 Theoretical Background 

In experimental measurements, errors may be divided into 
three classes: 

1. Blunders 

These errors are due to human mistakes, and all measure- 
ment programs suffer from them. Incorrect reading of meters, 
computer card punching mistakes, careless mating of waveguide 
flanges, etc., typify these blunders. In most cases, these 
errors will eventually be discovered and corrected. 

2 . Systematic Errors 


Systematic errors relate to the accurac y of a measurement 
and stem from incorrectly calibrated instruments, human meter 
reading error (e.g. , meter parallax), or such environmental 
conditions as strong magnetic fields interacting with meter move 
ments. Systematic errors will thus cause a measured value to be 
consistently high or low from the absolute value. Instrument 
calibration and the careful use of these instruments as recom- 
mended by the^manufacturer will minimize these errors. 

If a quantity is deduced from its relationship f = 
f(X^, X 2 , X 3 ,.... f X^) on n independently measured quantities 
(X^, X 2 , X^ X^) , then the worst case systematic error 
(accuracy) in f is 


n 

Af = l 
, i=l 


3f 

9X. 

l 


AX ± 


(5-6) 
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where AX . is the accuracy of the measured quantity X . . 
3 . Random Errors 


Random errors control the precision of a measurement or the 
number of significant figures that may be quoted. Such errors 
are noise-related and may stem from short-term gaussian noise or 
shot noise or they may be caused by the relatively long-term 
gam instability of an amplifier. If the variances of the 
quantities X^ (see above) are known, then the variance of the 
quantity f is given by the quadrature relationship 


2 

°f 


n 


l 

i=l 



5.5.2 PMIS Antenna Loss Errors 


(5-7) 


For PMIS, the waveguide and antenna are considered as a uni-! 
so that (3-2) reduces to 


L = 



(5-8) 


1 . Random Error 


It can be shown that the random fluctuations in the antenna 

temperature (a_ ) and in the sky temperature (o ) at X-Band are 
„ A ' L S 

two to three orders of magnitude smaller than the fluctuations in 

the uncorrected brightness temperature. Thus, from (5-7) , 


3L a 

3t b t b 


T - T 
A S 


T, 


A -“B -B (T a - T£)“ -B 
The uncorrected brightness temperature is computed by 


(5-9) 


where X = 




T ' = Tf + ( AT) X 


(5-10) 
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where T. and AT are constants established by laboratory calibration 
and X is the PCM count ratio. Thus, 


: “ t"?, + x Mt + (4T)2 %1 


1/2 


(5-11) 


The noise process responsible for a is short-term (within-the- 
minute fluctuations) and originates from the receiver front-end, 
from antenna phase shifters and from the residual PCM transients. 
It has been found that with a one-minute integration time, 


(AT) o 

x 


0.07 K (vertical) 

0.22 K (horizontal) 


(typically) (5-12) 


Fluctuations in T^ and AT originate from medium-term (within-the- 
week) gain instability as influenced by the synchronous demodulator 
in the receiver rear end. Repeated laboratory calibrations during 
the PMIS loss measurement period at A-Mountain established that 


' 1.60 

K 

(vertical) 

. 3.00 

K 

(horizontal) 

0.29 

K 

(vertical) 

1.05 

K 

(horizontal) 


(5-13) 


(5-14) 


as being typical. Using X = 3.50 (a maximum) and substituting 
5-12, 5-13 and 5-14 into 5-11, 


o 


T. 


B 


1.90 K (vertical) 

4.70 K (horizontal) 


(5-15) 
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Under typical conditions, 

T a = 288 K 
T s = 5 K 

(5-16 ) 

Tg = 120 K 
L a = 1.684 

Substituting 5-15 and 5-16 into 5-9, 


a 


0.019 (vertical) 
0.047 (horizontal) 


(5-17) 


The three-point convolution process used to borrow statistical 
information from neighboring beam positions decreases the noise 
on L a by 58%, giving a final random error 



0.010 

(vertical) 

°L = ' 


(total random error) 



0.027 

(horizontal) 


2 . Systematic Error 

r’-* 

The major source of systematic, or calibration error is the 
inaccuracy of the laboratory standard X-Band cooled load. Other 
errors, considerably smaller, include the accuracy of the sky 
temperature and of the antenna kinetic temperature. 

The worst case systematic error may be calculated by 


iL A - E 


T* - T 
B S 


<T a - 


2 * T A + 


t a - T S 

(T ft - T-) 


T 1 

A B 


1 


AT S 4- 


(5-19) 
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It is estimated that: 

AT = + 0.2 K (limited by accuracy of thermistors) 

ATg - + 0.4 K (limited by accuracy of knowledge in 

vapor volume absorption coefficient in 
accuracy of radiosonde instruments - 
see Section 4.4. 


AT' = + 1.0 K (See Section 5-2) 

£j — 


Substituting 5-16 and 5-20 into 5-19, 


aL A = ± (O * 0014 

sky temp, 
error 


0.0008 

ant. temp, 
error 


+ 0 . 0100 ) 

ref, cold 
load error 


(5-21) 


or 


AL a = + 0.012 


(total systematic error for vert, and 
horiz . ) 


3 . Total Uncertainty 

The total uncertainty is found by adding the random and 
systematic error, i.e. 


= i <0 LA + l A lJ ) 


(5-22) 


or 


E = 

+ 0.022 

(vertical) 

LA 

+ 0.039 

(horizontal) 


(5-23) 




i 
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5.5.3 PMIS Radome Loss Errors 

Since the waveguide loss is subsumed within the antenna 
loss for PMIS, (3-4) reduces to 

T - T 

T _ R S 

L R ~ + T, (L, - 1) - L„T’ ^ 


R A A 
1 . Random Error 


A B 


It can be demonstrated that random fluctuations in the 
brightness temperature have a negligible effect on (5-24) in 
comparison to the contribution from the noise on as computed 
from (5-18). Also, random fluctuations in the sky temperature 
and the radome temperature have a negligible effect on L . Thus, 
carrying out the indicated operations of (5-7) , using (5-24) , 


T - T 
R S 


l r (I A - T^)L A 2 l a 


(5-25) 


which, under typical conditions, reduces to 


o T = 0.6 a T 

l r l a 


(5-26) 


2. Systematic Error 


It may likewise be shown that the systematic error in L 
is, under typical conditions, given by 


R 


ALr . 0.6 AL a 


(5-27) 
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3 . Total Uncertainty 

It follows from (5-26) and (5-27) that the total uncertainty 
in the radome loss is 


e 


L R 


0.5 


which# using (5-23) , yields 


E 


l r 


+ 0.014 (vertical) 

+ 0.039 (horizontal) 


5.6 Effect of Errors in L^, L R on PMIS User 


(5-28) 


(5-29) 


What is the significance of these errors in the measured 
values from the viewpoint of the PMIS data user? In flight/ 
the antenna may be viewing a ground target or possibly a sky 
target? furthermore, it is covered by a radome so that (5-24) 
holds, except that T now represents the target (or scene) 
brightness temperature, and it is Tg that the user wishes to 
measure. Solving for Tg from (5-24) , 

T S - L A L R T 3 " V L A - l)T A - (L R - X > T R (5 - 30) 

5.6.1 Random Errors 


Even though there were random errors in the measured values 
of -and L R , the values themselves do not fluctuate. Thus, we 
may say that we know these values only to a certain precision. 
However, there is an appreciable random noise-caused error in 
T^ (more commonly known as the receiver sensitivity) , typically 
1 K over a short period of time. Long-term receiver instability 
may cause this value to become much higher over a period of 
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hours. Fluctuations in and T R are negligible. Applying {5-7) 
to (5-30) , the variance in the source brightness temperature is 


1 - [L R (T fl - T-nV + [I**B - <L a - 1)T a - T R l 2 a2 

o A K 




(5-31) 


B 


Assuming a source brightness temperature of 200 K, antenna and 

radome temperatures of 288 K and 293 K respectively and antenna 

and radome losses of 1.684 and 1.100, the uncorrected brightness 

temperature T' would be 241 K. The corresponding variances are 
a 


a 



a 



a 


2 

l r 


(0.010) 2 

vertical 

(0. 027) 2 

horizontal 

(l) 2 

vertical 

(l) 2 

horizontal 

(0 . 006) 2 

vertical 

(0.016) 2 

horizontal 


(5-32) 


Substituting these figures into (5-31) , 



0,94 + 0.25 

. 1.95 + 1.81 


+ 3.43 (vertical) 

(5-33) 

4- 3.43 (horizontal) 


antenna 

loss 

imprecision 


radome receiver 

loss noise 

imprecision 


or 


o 


2.2 K (vertical) 

2.7 K (horizontal) 


(5-34) 
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These values specify the precision with which a measured source 

brightness temperature may be quoted. They are larger than the 

receiver noise related system sensitivity because of the 

imprecision of the loss numbers. If the statistical fluctuations 

in T,-, are assumed to be normally distributed, then there is a 

probability of 68% that the measured value of T g will be within 

+ o of the mean. If the random error in and L R were reduced 

to zero, then cr_ = 1.9 K which is then the tangential sensitivity 

S 

of the radiometer system. 

5.6.2 Systematic Errors 

The systematic error, or accuracy, associated with T g is 
computed by using (5-6) and (5-30) in a worst case sense, yielding 

= IV T A - T B>I 4L a + l L A t_T A + T A> + i t a - V I al r 


A A B' 


'A R' 1 R 


+ W t b 


(5-35) 


The accuracies involved are 


1 Ah* 

— + 

0.012 

A 



1 AL-, 

= + 

0.007 

R 



j AT * 

= + 

1 K 

! B 



} so that, using 

the 

same 


(5-36) 


= + [0.62 


+ 0.59 


+ 1.85 


3 (5-37) 


antenna radome receiver 

loss loss calibration 

inaccuracy inaccuracy inaccuracy 
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or 


ATg = + 3 K (5-39) 

This means that the measured source brightness mean temperature 
Tg is within + 3 K of the true value. 

5.6.3 Total Uncertainty 

The total uncertainty in T g is found by adding the random 
and systematic errors, i.e.. 



or 


e 


+ 


5.2 K 
5.7 K 


(vertical) 

(horizontal) 


(5-41) 
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CHAPTER VI 


MFMR LOSS MEASUREMENTS 


6 . 0 INTRODUCTION 


The MFMR operates at L-Band (1.4135 GHz), K^-Band (18.0 GHz), 
K-Band (22.05 GHz) and K -Band (37.0 GHz) with fixed beam 

d 

antennas so that radiometric profiles in the flight direction are 

provided. The L-Band antenna is a linearly polarized stripline 

array, manufactured by AIL, with a 16° beamwidth and a beam 

efficiency of 95%. In order to obtain both vertically and 

horizontally components of incident radiation at L-Band, the 

entire antenna assembly is mechanically rotated through a roll 

angle of 90°. The K, K , and K-Band antennas are dual-polarized 

scalar horns, with simultaneous vertical and horizontal outputs 

to two radiometer channels for each band. 3 dB beamwidths are 

4.0° at K , 4.3° at K and 4.5° at K -Bands. The seven channels 
u a 

all use Hach [1968] radiometers with 100 ms integration times 
(in flight) and bandwidths of 27, 200, 200 and 500 MHz at L, K u , 

K and K -Bands respectively. Radiometric outputs and housekeeping 

cl 

data are all PCM encoded in a similar format to that used by 
PMIS. 


All antennas and radiometers are mounted on a positioning 
ring so that both roll (0° or 90°) or pitch (0° to 180°) motions 
are mechanically controlled. 

6 . 1 Mounting Configuration 

Fig. 2-18 shows the MFMR on the antenna positioner and tower 
inside the bucket. Luring the measurement sequence, the antenna 
positioner was set to compensate for MFMR pitch motions so that 
the beam was always pointing toward zenith. However, this caused 
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the position of the antennas relative to the bucket walls to 
change considerably, as shown in Fig. 2-7, with the antennas 
being very close to the top of the bucket when the pitch was 
90°. It will be shown later that this displacement had no 
measurable effect on the loss values. 

In flight, microwave absorber is placed on the aircraft 
bulkhead for two principal reasons: (1) the L-Band array is so 

close to the bulkhead for a pitch of 0° that mutual coupling 
between antenna and radome via bulkhead reflections must be 
eliminated, and (2) when the K u , K and K^-Band horns are pointing 
skyward and viewing a cold and constant sky, considerable radio- 
metric variations may occur in flight due to the backlobes viewing 
a variable terrain brightness temperature, unless this effect is 
artifically removed by causing the sidelobes and backlobes to view 
a relatively hot absorbing bulkhead. The absorber used is Emerson 
and Cuming flat Eccosorb, tuned for use at L-Band. 

In testing at A-Mtn, however, the Eccosorb was not furnished 
and PSL used Its own available absorber, the pyramidal material 
seen in Fig. 2-18. Unfortunately , this absorber is designed for 
use at X-Band and above and is essentially transparent at L-Band. 
The result is that all L-Band tests were with a non-absorbing 
bulkhead and thus did not simulate flight conditions. This 
situation will be discussed in detail in Sec. 6.5.4. 

6.2 Radiometer Receiver Calibration 


Each of the seven radiometer receivers was calibrated by the 
use of laboratory standard hot loads and cold loads. From (5-1) , 
letting T^ = ) and solving 2 simultaneous eqns. for the 

2 unknown calibration constants. 


T i = 


"A - T A 

X h - X C 


[ 6 - 1 ) 
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T, - T 

AT = ^- - X C (6-2) 

h c 

where X refers to the count ratio (see eqn. 5-2) and the sub- 
scripts h and c refer to hot and cold. 

During the testing phase, laboratory calibrations were done 
twice: (1) at the beginning of the MFMR tests on Feb. 8, 1975 

and (2) at the end of the MFMR tests on Feb. 21 , 1975. Table 6-1 
and Fig. 6-1 compare these constants over this period of time. 

No final calibration at K-Band was possible since the receiver 
failed during the intervening time. 



Table 6-1 

MFMR Laboratory Calibration 

Summary 



Initial Cal 

. (2-8-75) 


Final Cal. 

(2-21-75) 

Band 

T i < K > 

AT (K) 


T a. < K > 

AT (K) 

V 1 

329.87 

-54.24 


325.51* 

-54.48* 

V 2 

327.65 

-54.22 


323.11* 

-55.39* 

K a- 1 

330.57 

-58.73 


330.25 

-58.17 

K a“ 2 

331.14 

-60.02 


331.05 

-59.89 

K-l 

329.64 

-56.67 


** 

*•* 

K-2 

328.61 

-50.20 


** 

** 

L 

323.33 

-51.15 


323,84 

-50 . 93 


NOTES : 

* Final K u -Band calibration may be slightly inaccurate since 
there was insufficient equipment warmup time. 

** K-Band receiver was not operating at the time of final 
calibration-. — — - 


'V 
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6.3 MFMR Data Flow 


Loss values were measured at roll angles of 0° and 90° for 
the following pitch angles: 0°/ 5°, 10° , 15% 20% 25 °, 30°, 35°, 

‘0°, 45°, 50° , 60°, 70°, 80°, 90°, 100°, 110° r 120°, 130° r 140°, 
150% 160% 165% 170% 175° and 180°. In flight, a 0° pitch 
angle has all antenna beams directed toward nadir, with the L-Band 
array nearest the bulkhead and the horns nearest the radome nose. 

A 180° pitch angle corresponds to the zenith direction, with the 
horns nearest the bulkhead and the L-Band array nearest the nose. 

This arrangement thus requires 364 antenna loss measurements 
plus 364 more for each of the two radomes supplied, giving a total 
of 1092 measurements required. Most of these loss values were 
measured independently twice (and subsequently averaged) , so that 
roughly 2C00 measured values were processed. 

Due to an unforeseen interaction between the L-Band antenna 
and the radome, the radome loss values measured at that frequency 
are not valid, although an approximate value of L R = 1.09 
(independent of pitch angle) can be used. This is discussed in 
detail in Sec. 6.6. 

Antenna and radome loss values were calculated by the use 
of eqns. (3-2) and (3-4), with the sky brightness temperature 
(T ) estimated using the algorithm discussed in Chap. 4. The 
uncorrected brightness temperature Tp measured for each beam 
position was a 1-minute average. The 100 ms. integration time 
rms noise level (std. deviation) was monitored for each channel 
and a threshold, receiver sensitivity level of 6 K was established. 
This corresponds to a sensitivity of 0.25 K when referred to a 
1-minute integration time. Ti values noisier than this were 

JJ 

discarded . 


6.4 MFMR Antenna and Radome Loss Values (K._ , K, K -Bands) 


The 980 final loss values, as derived by this technique, are 
tabulated in Tables 6-2, 6-3 and 6-4. The antenna loss values 
are plotted vs. pitch angle in Figs. 6-2, 6-3 and 6-4. Cor- 
responding radome loss values for both radomes are plotted in 
Figs. 6-5 - 6-10. 

All of these loss measurements were made with PSL- furnished 
absorber in place on the bulkhead. The rise in L^ as the pitch 
angle approaches 180° (for a roll of 0°) is due to an increasing 
contribution to the apparent source brightness distribution T g 
as the horns come closer to the hot absorbing material. For a 
roll of 90°, this distance change with pitch angle is much less 
so that the effect is less pronounced. Radome losses are generally 
higher near 0° or 180° for a roll of 0° and are occasionally 
quite high from 170° - 180° since the antennas are looking 
straight up and into the radome fair ing-bu lithe ad interface, with 
the resulting likelihood of strong mutual coupling. 

It should be emphasised that the jagged nature of the loss 
graphs for L R is not due to errors in the measurement method, 
since these values were easily repeatable from one measurement 
set to the next under identical circumstances. Rather, this 
behavior stems from either resonant scattering between the horns 
and radome or resonant thickness effects (especially at K -Band) 

cl 

of the radome material. 

6.5 interaction Mechanisms between the MFMR Horns and the 


Aircraft Bulkhead 
6.5.1 Introduction 


As mentioned previously, it was observed that the antenna 
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Table 6-2 



PITCH 

L 

K -1 
u 

K -2 
u 

K-l 

0 ° 

1.365 

1.123 

1.119 

1.096 

5 ° 

1.366 

1.122 

1.119 

1 . 092 * 

10 ° 

1.365 

1.122 

1.118 

1.086 

15 ° 

1.364 

1.122 

1.118 

1 . 085 * 

20 ° 

1.363 

1.122 

1.117 

1.083 

25 ° 

1.363 

1.122 

1.117 

1 . 083 * 

30 ° 

1.361 

1.122 

1.117 

1.083 

35 ° 

1.361 

1.123 

1.118 

1 . 087 * 

40 ° 

1.362 

1.123 

1.118 

1 . 089 * 

45 ° 

1.361 

1.125 

1.119 

1 . 0B9 * 

50 ° 

1.362 

1.127 

1.121 

1 . 091 * 

60 ° 

1.360 

1.135 

1.129 

1 . 094 * 

70 ° 

1.360 

1.142 

1.134 

1 . 098 * 

80 ° 

1.361 

1.150 

1.143 

1 . 101 * 

90 ° 

1.360 

1.159 

1.153 

1.103 

100 ° 

1.360 

1.168 

1.162 

1 . 108 * 

110 ° 

1.360 

1.175 

1.168 

1 . 111 * 

120 ° 

1.360 

1.185 

1.175 

1 . 114 * 

130 ° 

1.360 

1.192 

1.181 

1 . 118 * 

140 ° 

1.360 

1.197 

1.187 

1 . 121 * 

150 ° 

1.361 

1.201 

1.193 

1.124 

160 ° 

1.362 

1.204 

1.196 

1.132 

165 ° 

1.362 

1.206 

1.201 

1 . 133 * 

170 ° 

1.363 

1.211 

1.204 

1 . 135 * 

175 ° 

1.365 

1.216 

1.208 

1.3 36 

180 ° 

1,366 

1.228 

1.219 

— 


1.082 1,124 

l . OBO * 1.123 
1.078 1.123 

1 . 077 * 1.121 

1.077 1.122 
1 . 078 * 1.122 

1.078 1.122 

1 . 081 * 1.122 
1 . 082 * 1.121 
1 . 084 * 1.122 
1 . 086 * 1.122 
1 . 089 * 1.125 
1 . 092 * 1.127 
1 . 096 * 1.131 
1.097 1.133 

1 . 102 * 1.137 
1 . 105 * 1.141 
1 . 109 * 1.144 
1 . 111 * 1.146 
1 . 115 * 1.149 

1.116 1.151 

1.124 1.156 

1 . 125 * 1.160 
1.126 1.164 


1.127 1.167 1.156 

1.173 1.160 



■*) 





ROLL = 

90 ° 



L 

rH 

> 

K -2 
u 

K-l 

K-2 

K . -1 
a 

K -2 
a 

1.369 

1.136 

1.140 

1.105 

1.098 

1.129 

1.114 

1.365 

1.135 

1.141 

1 . 104 * 

1 . 098 * 

1.128 

1.113 

1.363 

1.136 

1.139 

1.103 

1.098 

1.129 

1.113 

1.362 

1.136 

1.141 

1 . 103 * 

1 . 098 * 

1.129 

1.114 

1.362 

1.136 

1.142 

1.102 

1.098 

1.129 

1.113 

1.363 

1.137 

1.142 

1 , 102 * 

1 . 097 * 

1.130 

1.115 

1.362 

1.138 

1.143 

1.102 

1.096 

1.130 

1.115 

1.363 

1.140 

1.144 

1 . 100 * 

1 . 095 * 

1.130 

1.115 

1.363 

1,141 

1.144 

1 . 099 * 

1 . 094 * 

1.131 

1.116 

1.364 

1.141 

1,144 

1 . 098 * 

1 . 093 * 

1.131 

1.117 

1.363 

1.143 

1.146 

1 . 097 * 

1 . 092 * 

1.132 

1.117 

1.365 

1.145 

1.146 

1 . 096 * 

1 . 090 * 

1.133 

1.120 

1.365 

1.147 

1.148 

1 . 094 * 

1 . 089 * 

1.135 

1.121 

1.367 

1.146 

1.148 

1 . 092 * 

1 . 088 * 

1.136 

1.122 

1.367 

1.146 

1.148 

1.090 

1.087 

.137 

1.123 

1.369 

1.145 

1.147 

1 . 089 * 

1 . 085 * 

.137 

1.124 

1.367 

1.142 

1.145 

1 . 088 * 

1 . 084 * 

1.139 

1.125 

1.362 

1.138 

1.143 

1 , 087 * 

1 , 082 * 

1,138 

1.124 

1.365 

1.135 

1.141 

1 . 085 * 

1 . 080 * 

1.140 

1.126 

1.362 

1.133 

1.139 

1 . 083 * 

1 . 079 * 

1.141 

1.126 

1.365 

1.131 

1.138 

1.081 

1.077 

1.141 

1.127 

1.365 

1.131 

1.139 

1.082 

1.076 

1.147 

1.131 

1.366 

1.133 

1.141 

1 . 083 * 

1 . 078 * 

1.151 

1.135 

1.367 

1.134 

1.143 

1.084 

1.080 

1.156 

1.140 

1.372 

1.136 

1.144 

1.085 

1.079 

1.159 

• 1.144 

1.376 

1.138 

1.147 

— 

— 

1.163 

1.148 
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cm 


HFMR 

RADOME 1 LOSS (L„) 

K 






ROLL = 

0° 






ROLL = 

90° 



PITCH 

L 

K -1 
. u 

K -2 
u 

K-l 

K-2 

K 

a 

V 2 

L 

V -1 
u 

V 2 

K-l 

K-2 

V 1 

K -2 
a 



1.458 

1.300 

1.221 

1.217 

1.548 

1.521 


1.172 

1.161 

1.370 

1.293 

1.539 

1.574 



1.403 

1.276 

1.204 

1.198 

1,535 

1.512 


1.170 

1.160 

1.245 

1.212 

1.568 

1.591 



1.364 

1.252 

1.205 

1.190 

1.538 

1.531 


1.169 

1.160 

1.244 

1.213 

1.561 

1.578 

15° 


1-350 

1.242 

1.191 

1.176 

1.543 

1.549 


1.171 

1.159 

1.248 

1.223 

1.559 

1.555 

20° 


1.325 

1.234 

1.198 

1.183 

1,522 

1.536 


1.173 

1.161 

1.246 

1.220 

1.566 

1.545 

25° 


1.298 

1.223 

1.226 

1.209 

1.495 

1.504 


1.173 

1.162 

1.221 

1.199 

1.566 

1.536 

30° 


1.274 

1.216 

1.232 

1.212 

1.503 

1.527 


1.180 

1.165 

1.216 

1.193 

1.566 

1.530 

35° 


1.246 

1.208 

1.205 

1,185 

1.513 

1.542 


1.192 

1.170 

1.207 

1.187 

1.539 

1.512 

40° 


1.240 

1.206 

1.201* 

1.182* 

1.486 

1.512 


1.198 

1.176 

1.202* 

1.183* 

1.489 

1.463 

45° 


1.230 

1.194 

1.198* 

1.178* 

1.477 

1.50G 


1.212 

1.184 

1.197* 

1.178* 

1.483 

1.458 

50° 


1.222 

1.186 

1.194* 

1.175* 

1.472 

1.496 


1.205 

1.181 

1.192* 

1.174* 

1.474 

1.450 

60° 


1.214 

1.177 

1. 191* 

1.171* 

1.502 

1.518 


1.192 

1.174 

1.186* 

1.169* 

1.494 

1.481 

70° 


1.214 

1.193 

1.187* 

1.168* 

1.450 

1.466 


1.195 

1.186 

1.181* 

1.165* 

1.463 

1.440 

80° 


1.217 

1.209 

1.184* 

1.164* 

1.478 

1.487 


1.206 

1.204 

1.176* 

1.160* 

1.468 

1.454 

90° 


1.215 

1.207 

1.180 

1.161 

1.454 

1.461 


1.189 

1.175 

1.171 

1.156 

1.449 

1.436 

100° 


1.211 

1.204 

1.189* 

1.170* 

1.426 

1.429 


1.193 

1.177 

1.178* 

1.161* 

1.435 

1.425 

110° 


1.206 

1.199 

1.196* 

1.179* 

1.412 

1.414 


1.207 

1.188 

1.184* 

1.166* 

1.420 

1.413 

120° 


1.202 

1.194 

1.205* 

1.189* 

1.338 

1.340 


1.195 

1.176 

1.191* 

1.171* 

1.33B 

1.342 

130° 


1.208 

1.198 

1,214* 

1.198* 

1.354 

1.359 


1.171 

1.160 

1.198* 

1.175* 

1,331 

1.332 

140° 


1.214 

1.205 

1.223* 

1.207* 

1.357 

1.360 


1.163 

1.156 

1.204* 

1.180* 

1.345 

1.343 

150° 


1.207 

1.205 

1.234 

1.216 

1.362 

1.359 


1.165 

1.156 

1.211 

1.185 

1.367 

1.364 

160° 


1.292 

1.278 

1.342 

1.325 

1.440 

1.460 


1.177 

1.163 

1.271 

1.244 

1,401 

1.386 

165° 


1.378 

1.359 

1.420 

1.419 

1.486 

1.534 


1.190 

1.172 

1.383 

1.346 

1.409 

1.381 

170° 


1.379 

1.388 

1.446 

1.419 

1.552 

1.582 


1.195 

1.175 

1.417 

1.376 

1.445 

1.404 

175° 


1.412 

1.416 

1. 620 

1.595 

1.683 

1.770 


1.252 

1.209 

1.528 

1.463 

1.791 

1.698 

180° 


1.366 

1,323 

— 

— 

1.83B 

1.922 


6.880 

6.605 

— 

— 

9.473 

10.551 


cm 


*Interp elated Data 

NOTE : L-Band data missing because of mutual interact Ion effects between radotae and antenna 
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Table 6-4 

MFMR 

RADDME 2 LOSS (L„) 

K 


ROLL - 0° 


PITCH 

L 

V 1 

K ~2 
u 

R-l 

K-2 

K -1 
a 

K -2 
a 

0° 


1.536 

1.357 

1.241 

1.228 

1.520 

1.454 

5° 


1.475 

1.329 

1.228 

1.216 

1.577 

1.517 

10° 


1.426 

1.290 

1.231 

1.208 

1.601 

1.560 

15° 


1.408 

1.280 

1.230 

1.207 

1.595 

1.575 

20° 


1.394 

1.278 

1.232 

1.208 

1.571 

1.560 

25° 


1.367 

1.263 

1.219 

1.199 

1.544 

1.545 

30° 


1.339 

1.250 

1.212 

1.194 

1.547 

1.559 

35° 


1.301 

1.237 

1.210 

1.191 

1.531 

1.543 

40° 


1,291 

1.239 

1.214 

1.198 

1,569 

1.571 

45° 


1.278 

1.229 

1.189 

1.174 

1.515 

1.519 

50° 


1.229 

1.194 

1.153 

1.138 

1.482 

1.487 

60° 


1.229 

1.214 

1.157 

1.142 

1.459 

1.461 

70° 


1.235 

1.225 

1.162 

1.146 

1.512 

1.504 

CD 

°o 


1.239 

1.231 

1.171 

1.156 

1.463 

1.460 

90° 


1.248 

1.240 

1.205 

1.192 

1.457 

1.452 

J.00"' 


1.250 

1.243 

1.205 

1.199 

1.479 

1.468 

110° 


1.249 

1.245 

1.199 

1.183 

1.508 

1.523 

120° 


1.253 

1.246 

1.211 

1.190 

1.522 

1.527 

130° 


1.265 

1.256 

1.219 

1.195 

1.516 

1.518 

140° 


1.237 

1.230 

1.186 

1.173 

1.425 

1.447 

150° 


1.232 

1.235 

1.197 

1.189 

1.402 

1.419 

160° 


1.280 

1.279 

1.306 

1.299 

1.492 

1.546 

165° 


1.345 

1.358 

1.432 

1.435 

1.451 

1.490 

170° 


1.377 

1.404 

1.504 

1.511 

1.427 

1.489 

175° 


1.418 

1.469 

1.749 

1.766 

1.577 

1.664 

180° 


1.359 

1.422 

1.713 

1.738 

1.690 

1.774 


interpolated Data 

NOTE: L-Band data missing because of mutual interaction effects 
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ROLL = 

90° 




, K -1 

u 

K -2 
u 

K-l 

K-2 

K -1 
a 

V 2 

1.220 

1.196 

1.287 

1.224 

1.480 

1.537 

1.214 

1.193 

1.217 

■1,188 

1.463 

1.509 

1.207 

1.187 

1.228 

1.198 

1.478 

1.511 

1.194 

1.182 

1.239 

1,213 

1.495 

1.513 

1.185 

1.177 

1.237 

1.211 

1.494 

1.495 

1.180 

1.178 

1.237 

1.211 

1.487 

1-478 

1.180 

1.176 

1.236 

1.210 

1.486 

1.469 

1.185 

1.180 

1.229 

1.206 

1.486 

1.467 

1.191* 

1.185* 

1.223* 

1.202* 

1.481* 

1.463* 

1.197* 

1.190* 

1.218* 

1.199* 

1.476* 

1.460* 

1,203* 

1.195* 

1.212* 

1,195* 

1.471* 

1.456* 

1.20B* 

1.200* 

1.206* 

1.192* 

1.465* 

1.452* 

1.214* 

1.205* 

1.200* 

1.188* 

1.460* 

1.448* 

1.220* 

1.210* 

1.195* 

1,185* 

1.455* 

1.445* 

1.226 

1.215 

1.189 

1.181 

1.450 

1.441 

1.223* 

1.212* 

1.193* 

1.182* 

1.450* 

1.440* 

1.220* 

1.209* 

1.195* 

1.183* 

1.450* 

1.439* 

1.217* 

1.206* 

1.199* 

1.184* 

1.451* 

1.438* 

1.214* 

1.202* 

1.202*' 

1.184* 

1.451* 

1.437* 

1.211* 

1.199* 

1.206* 

1.185* 

1.452* 

1.436* 

1-208 

1.196 

1.212 

1.186 

1.452 

1.435 

1.222 

1.202 

1.243 

1.220 

1.436 

1.397 

1.228 

1.205 

1.287 

1.254 

1,393 

1.366 

1.238 

1.209 

1.294 

1.263 

1.468 

1.403 

1.282 

1.254 

1.411 

1.346 

1.885 

1.720 

1.319 

1.290 

1.313 

1,261 

1.924 

1.879 
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Fig. 6-2. MFMR antenna loss versus pitch angle — K u " Baild 
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Fig. 6-7. MFMR radome 1 loss 
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loss L_ at K-Band (and to a lesser extent, at K and K -Bands) 

rose monotonically with pitch angle increasing beyond 60° (when 

the roll was 0°) when the absorber was in place on the mockup 

bulkhead of the aircraft. When the absorber was removed, the 

antenna loss was essentially independent of pitch angle both 

at K and K -bands (no measurement at K-Band ware made without 
u a 

absorber). It is the intention of Sec. 6.5 tc show a mechanism 
of interaction between the bulkhead and the horn and to suggest 
further measurements which could improve the accuracy of the loss 
numbers applied to use in the operational situation on the P-3 
aircraft. Comments are restricted mostly to the K^-Band horn, 
but the approach can be generalized. 

6.5.2 Near-field Patterns of Circular Horns 


Assuming the scalar K^-Band horn has 25 dB sidelobes, the 
circular Taylor distributions for near-field given by Hansen [1964] 
for a 10 wavelength diameter aperture can be used by rescaling 
to the diameter of the K^-Band horn (36 cm) . Pig. 15 (p. 29) of 
Hansen lists several field patterns for a 10X diameter horn in 
free space. Let R be the distance to field point from the horn 
phase center (taken to be in the center of the mouth of the horn) . 

R 

Hansen defines A = — = — and plots field patterns for A = 0.037 5, 

2D Z /A 

0.05, 0.075, 0,125, 0.25 and «■. We are interested in the smallest 
two of these. Pig. 6-11 is a general plot of the power pattern 
(dB) for A - 0.0375 (R = 58 cm., 0.05 (79 cm.) and «> (far-field) 

TT Pi 

with the abscissa X = ka sin 0 = ~ sin 9 = 10ir sin 0. Since 
the K-Band horn has a diameter of 21. 7 X at 18 GHz, Pig. 6-11 

LI 

can be rescaled to this diameter, producing the graph of Fig. 6-12. 
According to Hansen, a typical circular Taylor design would require 
a monotonic aperture distribution with the power density at the 
edge of the horn iown about 8 dB from its value on axis. This 
produces the aperture distribution curve at a radius of 18 cm 
shown in Pig. 6-12. 
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Fig. 6-11. Fresnel and far-field patterns for 25 dB circular 
Taylor distribution. 
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The on-axis power density oscillates with distance from the 

phase center, assumed here to be at or very near the horn mouth. 

Fig. 6-13 is extrapolated from Hansen's data and plots the on- 

axis power in dB vs. distance. The variation out to about 200 cm 

is + 1 dB and can be seen to fluctuate most rapidly out to 50 cm 
~ 2 

from the horn. Beyond about 175 cm, the 1/R far- field behavior 
begins to predominate. 

6.5.3 Geometry of Horn Relative to Bulkhead 

Fig. 6-14 shows the geometrical relationship of the horn 
relative to the absorber-covered bulkhead, for pitch angles of 
180° (when horn is looking at zenith) and at 100° (slightly up 
from the aircraft flight vector) . The distance from the horn 
phase center varies from about 38 cm (closest approach at 0 = 180°) 
to about 81 cm (farthest approach at 6 = 70°) so that the near 
field pattern evaluated at the distance to the absorber must be a 
function of angle from horn axis as well as radius. To a first 
approximation, all near field patterns have the same on-axis gain 
level, and can therefore be normalized to 0 dB. 

6.5.4 Contribution of Absorber to Brightness Temperature 
The near-field curves of Fig. 6-12 make two assumptions: 

1) The horn is in free space 

2) Edge diffracted rays are negligible 

In reality, however, the horn is near a number of metallic surfaces, 
some of them very close. Furthermore, the approximate -8 dB 
edge taper introduces the possibility of strong near-field inter- 
action between the horn and the mounting ring, L-Band antenna and 
receiver enclosures. It has been shown in a series of papers by 
Peters, et. al. [1966, 1975] that a Fourier transformation of an 
aperture distribution which is assumed to be zero outside the 




Bulkhead 


Absorber 




Fig. 6-14. Positions relative to bulkhead - MF. 
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aperture of a horn (i„e., neglecting induced line sources at the 
horn edges) is satisfactory to predict the far- field shape of 
the main beam and the first few sidelobes, but fails to predict 
the level and structure of the far sidelobes and backlobes, 
either in the E-plane or in the H -plane. Typically, aperture 
integration methods which neglect the geometrically-diffracted 
rays from the edge give good agreement with measured results 
out to about 5 or 6 times the HPBW, when in the far field. 

Thus, in the near-field, a similar behavior can be expected 

so that the patterns of Pig. 6-12 may be considerably broader 

below about -20 dB, than shown. In addition, the in situ patterns 

may exhibit considerable asymmetry. 

v 

Can this near-field perturbed behavior, including edge 
diffracted rays, contribute to an increased antenna brightness 
temperature resulting from the 290 K absorbing wall? The answer 
is yes, and that it can be explained in terms of a broadened 
near-zone pattern behavior. The antenna brightness temperature 
is given by 

T a = JJ T B (0,<f»)G(0,<{,) dfl (6-3) 

4tt 

where T g is the source brightness temperature distribution 
(consisting of both the sky temperature and absorber temperature) , 
G is the antenna gain pattern, and dn = sin 6 d6 d<|> in a spherical 
coordinate system. This formula does not explicitly consider 
the distance from the antenna to the radiating surface. Normally, 
in the far-field, G(0,<|i) is independent of distance R from the 
antenna. However, in the near-field, the gain is a strong 
function of distance, i.e. 


G = G [6 , $ , R(0,<J>) J 


(6-4) 
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as shown, in Pigs. 6-11 and 6-12. A rough approximation to the 
K^-Band Brightness distribution (in the bucket) is shown in 
Fig. 6-15. Also shown are the far-field and estimated near- 
field patterns of the K u -Band horns, when the horn axis is at a 
pitch angle of 100° and a roll of 0°. In this situation, the 
phase center of the horn is approximately 71 cm from the nearest 
approach of the absorber. Edge diffraction from the horn and 
interaction with the adjacent K a and K-Band horns and the L-Band 
array may appreciably broaden and distort the near- field patterns 
from those shown in Fig. 6-12, leading to the dashed pattern 
(estimated) shown in Fig. 6-15 (top) . 

Denoting the pitch angle of the horn axis as 6 Q , eqn. (1) 
may be rewritten more explicitly as 

T A ( 0 O ) * 27 JJ T B <e\*')G(0 o - -4>) Sin e' de» d*' (6-5) 

4tt 

where the roll position is (<f» Q = 0°,90°) and (6' * $ 1 ) are the 
(pitch, roll) angles measured from the main beam position. This 
is recognized as the familiar convolution integral, and leads to 
an antenna temperature profile of the form shown in Fig. 6-15 
(bottom) . No scale units are placed on the ordinate since the 
function G(0 q - e f , <f> o - <f>') is not known exactly enough to carry 
out the integration in (6-5) . However, the gradual rise of the 
antenna temperature from 90° to 180° is due to the increasing 
contribution from the relatively hot absorber as e Q is increased. 

If the shape of the absorbing surface or its distance to the 
horn R(9 Q , <f> Q ) is changed by only a few cm, the near-field 
pattern of the horn can change appreciably so that the slope 

dT A 

in the near-zenith region (Fig. 6-15) may likewise change. 

°o 

Thus it is of considerable importance that the mockup aircraft 
bulkhead be as faithful a replica of the actual aircraft bulkhead 


gi 


PITCH ANGLE (DEGREES) 

Pig. 6-15. Source and apparent brightness temperature 
versus pitch angle. 
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as possible, especially when the ,? calibrated" antenna is to be 
used for near-zenith observation. 

When the absorber is removed, however, this critical relation- 
ship no longer exists and the antenna loss becomes nearly independ- 
ent of pitch angle, as shown in Fig. 6-16 for Channel 1 at K u ~Band. 
Similar behavior was observed at K -Band; no data without absorber 

c£ 

was available at K-Band since the receiver failed during that 
measurement sequence. 

6.5.5 Conclusions 

When the absorber is used on the aircraft bulkhead, the MFMR 
K^, K and K-Band loss values (both antenna and radome) measured 
at A-Mountain in February, 1975 are of questionable validity 
since the mockup and absorber used differ from the operational 
configuration. The complex nature of the interaction between 
the horn antennas, the sky temperature and the absorber temperature 
make it almost insuperably difficult to calculate the slope of 

the antenna loss and therefore the values of radome loss. 

o 

Future measurements should be taken with a much better bulkhead 
replica and should use the same absorber type and shape used 
in flight. 

6 . 6 MFMR Antenna and Radome Loss Values (L-Band) 


Table 6-1 lists the antenna loss values for the AIL L-Band 
antenna and shows that the loss (-1.34 dB) is essentially independ- 
ent of pitch or roll so long as there is no radome present or 
L-Band absorber on the bulkhead. 

The fact that the X-Band absorber has no effect at L-Band is 
illustrated in Fig. 6-17 for measured values, both with and 
without absorber. 
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It was found that this non-absorbing bulkhead produced a 
strong interaction between the L-Band antenna and the radome 
because of multiple internal reflections off the bulkhead. This 
effect is clearly evident in Fig, 6-18 where the uncorrected 
brightness temperature is plotted vs. pitch angle, both for 
radome (#2) on and off. When Tl (radome on) exceeds Tl 
(radome off) , the radome loss is greater than unity, but other- 
wise (0° < 0 Q < 5°; 90° <_ 0 Q < 165°) the radome loss would appear 
to be less than unity, an absurd result. Moreover, it is clear 
that the radome loss at L-Band should be essentially independent 
of pitch angle 6 o . The same effect is observed for the Aerojet 
L-Band array, as shown in Fig. 6-19, although the radome on 
temperature is consistently higher than the radome off temperature. 
This unstable behavior is again caused by internal scatter by the 
radome and is further evidenced in Fig. 6-20 where the VSWR of the 
Aerojet L-Band array is compared for both radome on and radome 
off. The peak mismatch ate 0 =5° agrees with the peak in T R at 

0 =5° in Fig. 6-19. 
o 


A partial correction of T ' due to the mismatch loss can be 
made even though it is not sufficient to describe the more general 
effect of redistributed antenna currents on the array. Thus, the 
radome loss would be estimated by 


T_ 

L R = T 


on 


off 


1 - 


2 ^ 


on 1 


1 ” i r off 


where T 


and T „„ refer to the curves of Fig. 6-19 
on off 

eted quantity is the mismatch loss correction. ]r] 
magnitude of the voltage reflection coefficient and 
to the VSWR curves of Fig. 6-20 by 


( 6 - 6 ) 

and the brack- 
is the 
is related 


_ VSWR ~ 1 
VSWR + 1 


(6-7) 


After L r is computed for all pitch angles 0 Q and averaged, a value 
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Fig. 6-19. MFMR Aerojet antenna temperature versus pitch 
anale — L-Band. 
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L r 8 1.09 (0.38 dB) (6-8) 

is obtained. It is suggested that this figure be used until new 
loss measurements (with a good L-Band absorber) become available. 

Since L-Band Eccosorb is being used in flight and must 
certainly contribute to the apparent source temperature, it is 
clear that this loss component was excluded from these measurements 
and that the values in Table 6-1 for L-Band are too low. 
Furthermore, (with absorber in place) would become even higher 
as 0 Q ^O 0 and the near-zone antenna pattern begins to intercept 
the absorber. 

6.7 MFMR Error Budget 

The general theory follows that used for PMIS, as discussed 
in Sec. 5.5.1. 

1. Random Error 

It can be shown that the standard deviations of the antenna 
and radome losses are given by 




T S )L W 


A 


[t a + V L w - 1} - 


l w t b^ 


'B 


and 


(6-9) 



[t r + V l a 


(t r - V 

1 > + t w l a (l w - 11 - Wb ]2 



(6-10) 


Under typical circumstances at K u , K, and K-Bands, T g = 13 K, 

T a = 288 K, T r - 290 K, T R * 42 K, = 0.25 K, L w = 1.05, 

L a = 1.065 and L R - 1.30. Using these values in (6-9) and (6-10), 


2. Systematic Error 


Systematic errors in the antenna loss arise primarily from 
inaccuracies in (1) the apparent sky temperature, AT s and (2) the 
receiver calibration, AT^. Systematic errors in the radome loss 
arise from not only AT and ATI but in addition, any systematic 
error AL^ in the antenna loss. The worst case systematic errors 
may be calculated from 


AL A T 


AT + AT’ 

A - T s S °A B 


{ 6-12 


and 


where 



* T s + l r l aV t b + ¥A 


(6-13 


d a = t a + T w (L w - - ¥b 

d r = t r + t a (l a " 1 > + T w L a [l w 1) Wb 

Using the assumed constants in 6. 7. 1-1, 

AL a = 0.004 AT„ + 0.004 AT' 

A S B 

and 


(6-14 

(6-15 

(6-16 


AL_ = 0.004 AT + 0.005 ATI + AL* 


(6-17 
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The accuracy of the receiver calibration, AT^, is related to 
the calibration of the reference hot and cold loads and again 
it is assumed that 

AT^ - + 1 K (6-18) 

The accuracy of the apparent sky temperature, AT„, depends not 
only on the adequacy of the Paris model for estimating the sky 
noise component, but also it depends very critically on the 
fidelity of the aircraft bulkhead mockup and its absorber geometry. 
Table 4-2 lists the systematic error in the sky noise component. 
However, there is apparently no way of estimating the systematic 
error due to the poor replication of the aircraft bulkhead 
environment in the mockup furnished. This becomes particularly 
critical at the pitch angles (0° - 20°, 160° - 180°) where good 
loss data is most needed. It is easily conceivable that a total 
systematic error AT = + 10 K could be encountered in such a 

i? * 

circumstance, thus yielding 

AL a = + 0.04 (6-19) 

AL_ = + 0.09 (6-20) 

XV * 

as a possible systematic error. 

This situation clearly leads to unacceptably high accuracy 
errors in the brightness temperature and can be improved only 
after a better replica of the aircraft bulkhead has been furnished. 
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CHAPTER VII 


BUCKET PERFORMANCE TESTS 


7.0 INTRODUCTION 

Several questions arise when evaluating the performance of 
the bucket and this chapter discusses the most obvious of these, 

e ,g. 


1. Does antenna location change within the bucket affect 
the observed brightness temperature? 

2. How does one know that the effective emissivity of the 
bucket doesn’t change over a period of many months or years? 

3. Does electromagnetic interference (EMI) affect operation 
in the bucket? 

7.1 Effect of Antenna Location Change 


It can be quite naturally suspected that the bucket, being 
a partially closed metal box, might have resonances associated 
with it. Clearly, as shown in Fig. 2-6, if the bucket is too 
small there are observable effects on even the "lit zone" portion 
of the antenna pattern. 

There is, however, clear evidence at L, K and K -Bands 
that location change within the bucket has a negligible effect. 
Figs. 7-1 and 7-2 compare the antenna loss vs. pitch angle at K^-2 
and K -1-Bands, with and without absorber on the bulkhead. Two 

a. 

conclusions are immediate: (1) the rise in with increasing 6^ 

(with absorber) is due to an increased noise contribution from the 
absorber (see also Sec. 6.5) and (2) when there is no absorber 
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ANTENNA LOSS L 



FITCH ANGLE (DEGREES) 

Fig. 7-2. MFMR antenna loss versus pitch angle — K a ~Band. 
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on the bulkhead, the flat loss curve indicates that raising the 
horns from a mid-bucket position (at 0 q =O° or 180° as in Fig. 2-7) 
to a position vary high in the bucket has a negligible effect. 

In early May, 1975 a similar experiment was conducted at 
L-Band using the S-194 radiometer receiver and a standard gain 
horn. The horn was connected to the receiver through about 
6 m of flexible coaxial cable so that it could be easily positioned 
at various locations within the bucket. The tests were conducted 
in mid-afternoon so that the sun was radiometrically visible by 
pointing the horn toward the' southwest corner of the bucket. 

Thus, some of the variation in T' was due to the sun intercepting 

a 

various portions of the antenna pattern. 

Even under these adverse conditions, however, a movement 
of the horn from near the bucket floor to near the top produced 
only a 1.7 K change. Pointing the horn directly toward the 
north wall of the bucket produced only a 3.9 K change. 

It may be concluded that if the radiometer has a reasonably 
low sidelobe level (say, below -25 dB) , there is virtually no 
dependence of the observed brightness temperature on antenna 
location within the bucket. 

7 . 2 Long-Term Bucket Emissivity Changes 

It has been pointed out in Sec. 2.4.2 that from a theoretical 
viewpoint even the most pessimistic of changes in the bucket 
emissivity or shape would not be observable by any modern radiometer. 

Nonetheless, certain baseline performance checks have been 
made at X, K, and K -Bands to produce an eventual verification 

tX a. 

that the bucket characteristics are stable over a period of 
months or years. The idea behind the verification check is that 
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if a calibrated radiometer in the bucket views a known source 
temperature Tg, then any change in the apparent loss would 
be due to a change in the emissivity or shape of the bucket. 

7.2.1 X-Band Measurements 


At X-Band, the sky temperature is typically 5 K and can be 
calculated (using radiosonde data) to an accuracy of +0.4 K. To 
ensure that the antenna characteristics would not change/ an 
X-Band standard gain horn (Scientif ic-Atlanta Model 21-8.4) was 
used and was mounted on the centerline of the bucket with the 
horn mouth 3/96 m above the floor, as shown in Fig. 7-3. The 
E-plane was oriented east-west. 

In order to bring the antenna brightness temperature within 
the normal operating range of the radiometer, an HP-Model X382 
precision attenuator was inserted between the horn and the radio- 
meter (the PMIS vertical channel was used) . The attenuator was 
set at 1.0, 1.2, 1.4 and 1.6 dB (four 1-minute averages each) 
and the corresponding values of T' noted (74.15, 82.78, 89.67, 

O 

and 97.87 K respectively). The sky temperature T was 4.9 K on 
the night of measurement (2-6-75), so that for the 1 dB setting, 

~ T s 296 - ^ 9 

l a = t“Tt? = m - -7- 4 - 1 * 312 (1 * i8 dB) (7 "^ 

A B 

which leaves 0.18 dB as the antenna loss of the X-Band horn. 
Continuing in this fashion, the graph of Fig. 7-4 is obtained, 
in which the uncorrected brightness temperature T^ is plotted 
vs. the total loss in dB. The measured points give a slope of 
4.0 K/0.1 dB and the calculated slope (obtained by adding 0.18 dB 
to the attenuator setting and using T sky = 4.9 K) gives 
4.8 K/0.1 dB. The difference in slope may be attributed to using 
the calibration constants T^, AT outside of their 100 K - 200 K 
intended range. 
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If the emrssivity of the bucket were to degrade, then T g - k T sk ^ 


where k (>1) is the degradation factor. The uncorrected 
brightness temperature is then given (from 2-7) by 


T 1 = T + 
A 


kT , - T,. 

sky A 


(7-2) 


where all quantities except k are now known. It is straightforward 
to show that a change of (by attenuation setting) from 1.5 dB 


to 1.6 dB causes a change in AT ' given by 

D 


T„ - kT 


AT . = 


sky 


(.033) 


(7-3) 


A 


where is the median value, = 1.429. Therefore, 


AT 


B 


(7-4) 


For example, when T ft - 296 K (thermistor reading) and T g - = 4.9 K, 
and assuming k = 1, 


ATg - 4.7 K/0.1 dB (7-5) 

If, due to bucket corrosion, deformations, etc. k were to double, 
then AT I = 4.6 K/0.1 dB. Extending this to a 1.0 dB change in 
L^, a 1 K difference in slope would be observed. 

Clearly, a change in bucket emissivity is indistinguishable 
from a systematic error in T^y* Therefore, it is imperative in 
constructing a history of the bucket performance that be 

accurately assessed from sounding data. 
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7. 2., 2 L, K u , K a ~Band Measurements 

Similar bucket performance tests were conducted at L, K u and 
K -Bands except that a 183 cm diameter microwave absorbing disk 
was suspended 6 m above the MFMR antenna assembly as shown in 
Fig. 7-5. The purpose of the disk was to present a known hot 
brightness temperature which subtended most of the antenna beam 
solid angle [Kraus and Carver, 1974] and thereby reduce the 
systematic error in the apparent sky brightness temperature. The 
disk was held rigidly in place by a quadrupod arrangement of 
four aluminum spars and the absorber temperature was averaged 
from the readings of eight thermistors embedded at various 
locations across the disk. 


Five one-minute averages were recorded for each of the 
five channels, with a pitch angle of 180° in all cases. The 
roll angle was set to 0°, then to 180° and then back to 0°, 
with the results as shown in Fig. 7-6. The sky temperatures 
(2-18-75) at L, K u and K a -Bands were 4.3 K, 6.3 K and 11.1 K 
respectively. The average disk temperature was 288.5 K. It is 
clear that the high T' values indicate that the main beam of 

•D 

each antenna was viewing primarily the disk and that errors in 



would be of negligible importance in determining the 


repeatability of the experiment. 


However, the lack of repeatability of the T^ in the two 
measurements for a roll of 0° indicates that the mechanical 
repositioning accuracy of the antennas was poor. The roll angle 
settings were carried out using hand positioning of the ring 
(by LEC/HASD personnel) and the lack of repeatability may have 
been due to personnel fatigue and consequent carelessness or 
it may have been an inherent problem in the MFMR positioning 
system. Until this problem is corrected, it is not reasonable 
to expect that these results can be repeated. 


ORIGINAL PAGE IS 
OP POOR QUALITY 



Fig 


7-5 


Absorbing disk located above MFMR antenna assembly 
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FREQUENCY 

PITCH = 

180°, ROLL = 0° 

PITCH = 180°, ROLL = 90° 


NO. 1 

NO. 2 


k u .gh. 1 

174.9K 

168.8K 

221 .9K 

K u -CH.2 

176.6K 

170.2K 

223.3K 

Kg - CH. 1 

191. IK 

184.9K 

197.9K 

K a -CH. 2 

1S3.8K 

187 .2K 

197.9K 

L 

188.1 K 

187.4K 

181.9K 


Pig. 7-6. L, K u , K a bucket verification tests. 
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7.3 EMI 


Systematic spectrum surveillance tests were conducted 
throughout the VHF-UHF bands to identify and record interfering 
RF signals that could possibly be detected by the radiometer 
IF circuitry. These tests made use of the HP 8550 series 
spectrum analyzer with various Yagi and log periodic antennas. 
The antennas were placed in and near the bucket and were pointed 
in different directions. 

The only interfering signal that seemed to be correlated 
with the radiometer output was a 150 MHz sporadically actuated 
carrier from a state police repeater tower located approximately 
600 m northeast of the bucket. This signal level occasionally 
exceeded an input power level of -20 dBm when a moderately 
directive three-element Yagi was pointed directly toward the 
repeater tower. 

No power line interference or ignition noise was observable 
at the radiometer output. 


APPENDIX 


PHIS & MFMR RADIOMETER PROGRAMS 

INTRODUCTION 


Two separate sets of program procedures are given here: one for the PUIS 

series of programs, and one for the MFMR series of programs. Although some 
repetition is entailed, the procedures for both series are written as separate 
modules, designed to be used independently of each other without cross-refer- 
encing. 

To facilitate use of the procedures, a reference index is provided for 
each series of programs. The pagination style is a key to the series of pro- 
grams (PMIS or MFMR), the program within the series (PMIS 1, PMIS 2, SKYTEKP 
program used with PMIS data, etc.). For example: Pl.l denotes page 1 of the 

PMIS 1 program procedures; PS.l is page \ of the SKYTEMP program as used in 
the PMIS series; and, P3.1 is page 1 of PMIS 3 program procedures. In like 
fashion, the prefix rl M n indicates a page of the MFMR procedures writeup. 
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PHIS PROGRAMS - REFERENCE INDEX 

Page No . * s 

GENERAL V.l - P.5 


A. To load and execute System-7 PCM Program ••>••••••<■ P.l 

B. Instructions of A., above, for newcomer to System-7 P.la 

C. Telemetry Playback console - sketch and power settings . ... P.2 

D. Comments on System-7 Program P.3 

E. Example: printout from control typewriter, PMIS data run, . . P.4 

F. Notes on control typewriter printout 

1. Explanation of terms P.5 

2. "Overflow" P.4,5 

3. "Valid-Data-Flag Drop". ......... P.4,5 

PMIS 1 PROGRAM Pl.l - PI. 3 


A. Card setup for running PMIS 1 on computer Pl.l 

B. Discussion of PMIS 1 card setup 

1. Phase 1 (PMIS 1) PI. 2 

2. Phase 2 (PMIS 2) PI. 2 

3. Phase 3 (Job-Step Disk-Dump), .............. PI. 2 

C. Optional Control Cards - PMIS 1 (1st phase) 

1. UPSI card PI. 2 

2. Disk Address card (as used in PMIS 1) PI. 3 

D. Optional Control Cards - PMIS 2 (2nd phase) 

1. // EXEC EM03B card; Tl and AT card PI. 3 

PMIS 2 PROGRAM P2.1 - P2.7 


A. PMIS 2 Program as 2nd phase of a PMIS 1 job setup 

1. "PMIS 2 Card" output by PMIS 2, input to PMIS 3 .... . P2.1 

2. Job-Step Disk-Dump P2.2 

B. PMIS 2 Program as single-phase, separate job 

1. PMIS 2: separate job or 2nd phase treatment. ...... P2.3 

2. Control cards: when optional, when mandatory ...... P2.3 

C. Card setup for running PMIS 2 (separate job) on computer. . . P2.4 

D. Control cards 

1. Tl and AT card P2.5 

2. Disk Address card (as used in PMIS 2) P2.5 

E. Notes on conversion to engineering units in PMIS 2 Program 

L. Housekeeping data .......... P2.6 

2. Uncorrected sky brightness temperature, T T g P2.6 

3. Standard deviation of uncorxected sky brightness 

temperature, o^ . P2.6,7 
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PHIS PROGRAMS - REFERENCE INDEX, «mfd 


Page No.’s 

SKYTEMP PROGRAM PS.l ~ PS. 3 


A. Input data cards for SKYTEMP program 

1. Header card, description and format PS.l 

2. Pressure, temperature, dew point (PTD) cards - format . . PS.l 

B. Card setup for running SKYTEMP on computer PS. 2 

C. "SKYTEMP Card" output by SKYTEMP, input to PMIS 3 PS. 3 

PMIS 3 PROGRAM P3.1 - P3.7 


A. Antenna loss factor, radome loss factor; computational 

equations P3.1 

B. Standard deviations: antenna loss, radome loss • P3.2 

C. Card setup for running PMIS 3 "Radome Is Off" on computer . . P3.3 

D. Discussion: card setup for "Radome Is Off" job deck 

1. Stacking data decks P3.4 

E. Control cards for "Radome Is Off" job deck 

1. "Radome Is Off" card. P3.4 

2. "PMIS 2 Card" P3.4 

3. "SKYTEMP Cards" P3.5 

F. Card setup for running PMIS 3 "Radome Is On" job on 

computer P3.6 

G. Discussion: card setup for "Radome Is On" job deck 

1. Stacking data decks ........... P3.7 

H. Control cards: "Radome Is On" card P3.7 



PROGRAM PROCEDURE TO REDUCE PHIS DATA TAPES 


PRIOR TO LOADING THE SYSTEM-7 PCM PROGRAMS 

PERFORM CHECKOUT OF ALL PCM HARDWARE AND TAPE(S). 


TO LOAD AND EXECUTE THE SYSTEM-7 PROGRAM, PROCEED AS FOLLOWS: 

(1) RUN U-ZERO TAPE (LOAD UNIT ADDRESS = X'0000') 

(2) LOAD IPL PROGRAM FROM DISK (LOAD UNIT ADDRESS « X'00Q2') 
*** I POO A ENTER CONTROL STATEMENT 

(3) CYCLE HOST ATTACH SWITCH TO (ENABLE AND IPL) AND 

THEN TO (ENABLE) j 

IT MUST REMAIN ON (ENABLE) FOR REMAINDER OF PROGRAM EXE- 
CUTION, 

(4) SET DISK PROGRAM "DATA SET NAME" BY FOLLOWING R(EFER) 
STATEMENT: 

R J0BLIB,F2,,PCMUSER 

(5) L(OAD) PCM PROGRAM INTO CORE BY FOLLOWING L(OAD) STATEMENT: 
L PCM01B 


THE SYSTEM-7 IS NOW PREPARED TO ACCEPT LINKUP WITH THE SYSTEM-370 
AND WILL REMAIN IN A WAIT STATE UNTIL PMIS1 PROGRAM IS EXECUTED 
BY THE SYSTEM-370. 
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PROGRAM PROCEDURE TO REDUCE PMIS DATA TAPES 


Notes: "PMIS" Is the name of the program and of the radiometer being tested; procedures 
are given in large type, while explanatory comments are given in smaller type enclosed 
by hand-drawn brackets, thus: Oor benefit of one who is unfamiliar with the System- 

PRIOR TO LOADING THE SYSTEM-7 PCM PROGRAM: 

PERFORM CHECKOUT OF ALL PCM HARDWARE AND TAPE(S) . 

TO LOAD AND EXECUTE THE SYSTEM-7 PROGRAM, PROCEED AS FOLLOWS: 

(1) RUN U-ZERO TAPE (LOAD UNIT ADDRESS = X'0000’) 

£u-ZERO, the tape to IPL the computer. Is a small, blue punch tape kept 
right by the cu„puter; put Address switches on 0000 (Addresses are con- 
trol switches on front of computer; there are 4 Address switches, thus 
each Address switch is set on 0, as Indicated by the four zeros in 
step (1), above 

(2) LOAD IPL PROGRAM FROM DISK (LOAD UNIT ADDRESS = X'0002') 

£S?L denotes "Initial Program Load". Leave the first 3 Address switches 
on zero, move the right hand Address switch to setting of 2.^| 

*** IPOO A ENTER CONTROL STATEMENT 

^pie above is a statement that the computer sends back to you after 
you have done step (2) 

(3) CYCLE HOST ATTACH SWITCH TO (ENABLE AND IPL) AND THEN TO (ENABLE) 
WHERE IT MUST REMAIN FOR THE REMAINDER OF PROGRAM EXECUTION. 

|j:he Host attach switch is another control switch on the front of the 
computer; switch it to "Off" (that’s down), then up to "Enable and IPL" 
and back to center position, which is "Enable"; it must stay on "Enable" 
all the time you're transferring data from Sys.-7 to Sys.-370 - this is 
done in order to establish a linkage with the System-370 .T 

(A) SET DISK PROGRAM "DATA SET NAME" BY FOLLOWING R(EFER) STATEMENT: 

j]|tep (4) tells you to type an "R", that is, "Reference Statement", on 
the control typewriter. The statement to be typed, is as follows, be- 
low (be sure to leave a space between the "R" and the "J")Q 

R J0BLIB,F2,,PCMUSER 

(5) L(OAD) PCM PROGRAM INTO CORE BY FOLLOWING L(OAD) STATEMENT: 

Jafter typing the statement of step (4)*s instruction, execute step (3) 
by typing the following statement, below. Leave a space between "L" 
and "P".3 

L PCM01B 


THE SY3TEM“7 IS NOW PREPARED TO ACCEPT LINKUP WITH THE SYSTEM-370 AND 
WILL REMAIN IN A WAIT STATE UNTIL PMISl PROGRAM IS EXECUTED BY THE 

system-370. 
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COMMENTS ON SYSTEM-7 PROGRAM 


PAGE 



3.TMT 


SOURCE STATEMENT DOS ASM/7 

(360A-TX-011) VI M2 

12/10/74 


2 

* 

PC MO 10 

PROGRAM SPECIFICATIONS! 


000000? 0 


•3 

* 

PURPOSE 

; INPUT (MFMR ) OR (PM IS) DATA INTO A 648 WORD SUFFERS 

00000040 


4 

* 


I A OR Bi, DETERMINE TYPE { MFMR J OR (P.MISl, CHECK FOR 

00000050 


5 

ft 


PROPER SEQUENCE OF (TIME-WORD- 

3— FLAG) AND (FIO) FORMS, 

00000060 


6 

3* 


TRANSFER TOGGLED BUFFER AT * TIME- WORD-3 ) IF (VALID DATA! 

00000070 


7 

ft 


FLAG IS PRESENT, AND SIGNAL ( ENO-OF-F IL El WHEN 

oooooooo 


8 

ft 


( VAL I D-DAT A— F L AG) DROPS. 


00000090 


9 

ft 




OOOOOiOO 


10 

ft 




OOOCOliO 


11 

ft 

BUFFER 

FORMAT J 


.00000 L 20 


12 

5* 



* 

0C000130 


12 

ft 


HEADER (8-WORDS) DISP (0000) 

CODE WORD (0000) DATA 

00000140 


14 

ft 



( FFFF) EOF 

00000150 


15 

ft 



(EEEE) EOJ 

00000160 


16 

ft 



BLANK 

00000170 


17 

ft 


(0002) 

LOGID 2-CHARS (EBCDIC) 

00000 Uio 


18 

ft 


(0003) 

LOGID 2-CHARS ( EBCDIC) 

00000190 


19 

ft 


(0004) 

TIME-WORDI 

00000200 


20 

ft 


(0005) 

TIMF-WQRD2 

000002 LO 


21 

ft 


(0006) 

TIME-WORD? (QUEUE XFR ) 

00000220 


22 

ft 


(0007) 

BLANK 

00000230 


23 

** 




00000240 


?4 

ft 


FRAME NUMBER U1 (0008) 

SYNC 1 

00000250 


25 

ft 


(0009) 

SYNC 2 

00000260 


26 

ft 


DATA11 (30 WORDS) (0010) 

DATA 

00000270 


27 

* 


TH D U 


0G000280 


28 

ft 


(0039) 

DATA 

00000290 


29 

ft 




00000300 


20 

ft 


FID! (0040) 

DI01 (FID CHK , SET TYPE) 

00000310 


31 

ft 




00000320 


22 

ft 


DATA12 (127 WORDS) (0041) 

DATA 

00000320 


33 

ft 


THRU 


00000340 


34 

ft 


(0167) 

DATA 

00000250 

M3 

25 

ft 




00000360 

(U 

CM 

ro 

36 

ft 


FRAME NUMBER (2) (0168) 

THRU (0169) SYNC (TIME CHK) 

00000370 


hi 

03 


EXAMPLE OF PRINTOUT FROM CONTROL TYPEWRITER - PHIS DATA RUN 


The following System-7 messages will he executed by the System-7 when the Sys- 
tem-370/System-7 linkup occurs, as indicated by the statement "S370 is ready 
to accept data". Responses typed by the operator on the control typewriter 
are indicated by enclosure in a box. 


*** IPOO A ENTER CONTROL STATEMENT 
j R JQBLIB,F2, ,PCMUSER ] 

L ECMOlin 

S370 IS READY TO ACCEPT DATA 

ENTER OPTIONS BY TYPEWRITER (REQUEST) : 

REQ (LOG) TO ENTER LOG-ID WORD THEN (OSCj OR (RUN). 

REQ (OSC) TO RUN SIMULATED PCM DATA. 

REQ (RUN) TO RUN PCM DATA. 

REQ (EOF) MARK FILE AND IGNORE DATA. THIS IS NORMALLY 
AUTOMATIC BY VALID-DATA-FLAG DROP. 

REQ (EOJ) TO SIGNAL TERMINATE JOB. 

NOTE — NEW FILES MAY BE OPENED BY (LOG) /(OSC) OR (RUN) REQUESTS 
AFTER AN (EOF) REQUEST. 


LOG 


ENTER 4-DIGIT (HEX) LOGBOOK(ID) WORD 
1 P019 


REQUEST (OSC) OR (RUN) TO PREP FOR PCM DATA 


RUN 


ENTER PCM DATA 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 


EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 


OR1: BUN 


ENTER PCM DATA 

EOF PROCEDURES EXECUTED 

EOF PROCEDURES EXECUTED 


{ see following page, "Notei on 

2£?s\ T r"s lter Frintout ■ 

PMIS Data Run 

job was started again at selected place on tape 


QRl; rEOJ~] -4 — "^2 — - — -z — end of last tape; end of job was requested 

S7/S370 LINK TERMINATED 
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NOTES ON CONTROL TYPEWRITER PRINTOUT - PHIS DATA RUN 


Explanation of Terms 

LOG ID = name of the tape (4-digit word). The first digit must be 
either "P" or "M" in order for the computer to recognize 
data as being either PMIS or MFMR data. 

OSC = running oscillator to de-bug tape. 

PCM = pulse code modulation. 

VDC = Test Id. number (name of file) . The Test Id. is the 
"valid data" code, or word. 

EOF - end of file. 

EOJ = end of job. 


Overflow 


When overflow occurs, as it did during the processing of PMIS tape 
F019 (indicated by arrow on "Example of Printout from Control Typewriter - 
PMIS Data Run" on the previous page) , determination of where to re-start 
the job is made by: counting EQF's; and, by checking the time on the time 
code reader and the time code on the log. The tape is then backed up to 
the selected point, and the job is started again. 

Although whole files are not lost prior to overflow ( all files are lost 
afterward), it is possible to lose 25 data cycles, spread throughout the 
files, out of the job. There are 640 words per data cycle. 

When data overflows: a bell rings; EOF messages cease to print out; 

and, the control typewriter skips lines where the messages would have been. 


Valid-Data- Flag Drop 

The twenty-second word of each PCM data cy ?a is the "valid data word" 
(VDC). The high-order bit of this word is the ’ Valid-Data-Flag". The 
valid-data-f lag is one (1) during each file and is zero (0) between files. 

The System-7 program monitors the valid-data-f lag and executes an end- 
of-flle procedure each time the valid-data- flag drops, i.e. , goes from 1 to 0. 
The PMIS I and MFMR 1 programs also monitor the value of the entire valid- 
data-word; in case the System- 7 fails to execute an end-of-file procedure 
any time the VDC changes, the PMIS I and MFMR 1 programs will do so. 


page P.5 



TAPE NAME & DATE CREATED 


TAPE LOC. 


DISK 


VOLSER NO. 


bRZOOl 


282 


PUNCH 


2701 SEL 


REAL ALLOCATION 


COMMENTS 


FORMS 1 2 4 


e - 7 Trans Pe i 
fcun Tn f-nre: ^LrOKH 


SPECIAL 


CARRIAGE TAPE 


950 Rev. 8-74 


page Pl.l 


ORIGINAL PAGE IS 
OF POOR QUALITY 


TAPE DR 


H 1 8 (PM1JS & /armrs/n c^arcL — e>tft te/ta' in 

(P'MZS 2, t - oft ten a/ //I PMI5 1 ) /* 


PMIS 1 PROGRAM CARD SETUP 


The PMIS 1 program is executed by running the following setup on the 
System-370: 
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v 7 joi: ' 


Priis 1 & 
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APPROX RUN TIME 
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PHIS 1 PROGRAM CARD SETUP & PROGRAM OPTIONS 


PMIS 1 Program Card Setup 

1. Phase 1 

The example of a "PMIS 1 Program Card Setup" given on the previous 
page actually includes three separate programs, PMIS 1, PMIS 2, and Job- 
Step Disk-Dump, all executed by the PMIS 1 program setup. When the three 
programs are thus used, PMIS 1, PMIS 2, and Job— Step Disk— Dump are the 
first, second, and third phases, respectively, of a single job. 

By removing the PMIS 2 cards numbered 12, 13, and 14, PMIS 1 can 
be run as a separate job. 

2. Phase 2 

PMIS 2 may also be run singly, by choice. However, If certain val- 
ues mus: be changed in PMIS 2 data run previously, there is no option; 
PMIS 2 must be executed as a separate job (see PMIS 2 writeup and card 
setup) . 

3. Job-Step Disk-Dump 

When PMIS 2 (or PMIS 1, if PMIS 2 was not included in the job setup) 
has processed all of the data on the disk, control is turned over to Job- 
Step Disk-Dump (cards #15 and #16). Job-Step Disk-Dump prints all of 
the data cycles that were dumped into the disk by PMIS 1. 

This job step also computes the amount of run time for the entire 
program and prints this time on the SYSLOG and on the printer, thus must 
be Included in all foreground PMIS and MFMR programs. Disk Dump dumps 
the disk only if it is enabled by following the If EXEC EM03H card (#15) 
with the dump card (#16). If no dump is desired, follow the // EXEC 
EM03H with the /* card (#17). 


Optional Control Cards - PMIS 1 (First Phase) 

1. Card No. 8, the DPSI card, is optional. It contains a string of 

eight binary digits. 

The binary value of the seven low-order bits of this string Is the 
number of data cycles which will be dumped onto the disk each time a new 
Test Id. is established, i.e., the number of data cycles dumped per file. 
Note: the last data cycle of each file (Test) is not under the control 

of the // DPSI card and is always placed (dumped) on the disk. 

The high-order bit is a flag to determine whether or not the diag- 
nostic is to be printed from PMIS 1 (Phase 1 of the example, "PMIS 1 
Program Card Setup"). If the high-order UPSI bit Is 1, the full print- 
out is obtained; if the high-order bit is aero, or if there is no UPSI 
card. Phase 1 printout is suppressed. 
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PHIS 1 PROGRAM CARD SETUP £ PROGRAM OPTIONS 


Optional Control Cards - FMIS (First Phase) , continued 

2. Card 10, "Disk Address" card, must be used with great caution. 

Its misuse can destroy all data on the disk, since it updates, i.e, , 
changes the disk directory. On the card shown in the PMIS 1 program 
setup, each word has the value "99", which is the starting address. 
This initializes the disk and causes it to start loading data at the 
beginning of the disk file. 

Format of Disk Address Card (when used in PHIS 1) : 3110. 

a. First word, HIGH(j)N: the highest disk address that PMIS 1 or 

MFMR 1 has operated on up to a given instant; 

b. Second word, HIGH(j>FF: the highest disk address that PMIS 2 

or MFMR 2 has operated on up to the same given instant; 

c. Third word, J(j)BEND: the highest disk address that the last 

PMIS 1 or MFMR 1 job has operated on up to the time it was 
terminated. 

Note: the disk control is set up such that both PMIS 1 and 

PMIS 2 can operate simultaneously in different partitions of the 
computer. And, if PMIS 1 is in process of executing while PMIS 
2 is also in the process of executing, then PMIS 2 has to know 
exactly where PMIS 1 is. Three variables are transmitted through 
the disk so that these two programs can communicate with one 
another; these variables are HIGH(j)N , HIGH<j)FF, and J(j)BEND. 


Optional Control Cards - PMIS 2 (Second Phase) 

3. Cards numbered 12 and 13, the "// EXEC EM03B" and the "T1 and AT" 

cards, respectively, are optional in the PMIS 1 program. Inclusion 
of these cards causes the PMIS 2 program to be executed in the PMIS 
1 program (see example of PMIS 1 program card setup). However, the 
T1 and AT card is optional, and it can be included only if the 
values of T1 and AT are known prior to execution of the job. If 
omitted, the values used will be those that were used the last time 
the PMIS 2 program was executed. 

Format of T1 and AT Data Card: 4F10.4. Words are: T1 vertical; 

AT vertical; T1 horizontal; and, AT horizontal. respectively. 

If it should later become necessary to change tin values of T1 
and AT in data that were run previously, this is done by re-running 
PMIS 2 as a single phase (s^.e PMIS 2 writeup). 
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PHIS 2 PROGRAM 


When the EOJ instruction is issued at the System-7 terminal, the link- 
age between the System-370 ana the System-7 is terminated, and program PMIS 1 
(the first phase of. the job) ends and turns control over to the second phase 
(PMIS 2 program) . 

PMIS 2 retrieves the raw data placed on disk DRZ3Q1 by PMIS 1 and com- 
putes the engineering units; it then places these results on the disk, prints 
them out, and punches a card for each file (Test). These cards are referred 
to as "PMIS 2 cards” (see example) . 

The "PMIS 2 card" will be used in the PMIS 3 program (see example of 
PMIS 3 program card setup) . 



Note ; the Log Id., above, is the name of the tape. The first digit must be 
either "P" or "M" in order for the computer to recognize data as being either 
PMIS or MFMR data. The Valid Data Word (or Valid Data Code, VPC) is the Test 
Id. number (name of file). 
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PHIS 2 PROGRAM 


When PMIS 2 has processed all of the data on the disk, control is turned 
over to Job-Step Disk-Dump. This is the last job step of the three-phase job 
illustrated by the PHIS 1 Program Card Setup. Job-Step Disk-Dump prints all 
of the data cycles that were dumped into the disk by PMIS 1. 


Note; the last data cycle of each file (Test) is not under the control 
of the // UPSX card and is always placed on the disk. 


Job-Step Disk-Dump is described in "PMIS 1 Program Card Setup & Program 
Options" under "PHIS 1 Program Card Setup". 


page P2. 


PHIS 2 PROGRAM - USAGE AS SINGLE-PHASE JOB 


Whenever it becomes necessary to re-run PMIS 2 program at a later time 
with new values of Tl and AT, PMIS 2 must be executed as a single-phase job. 

The PMIS 2 program card setup may include or omit optional control cards 
as follows: 

1. as Phase 2 . .of a three-phase job setup . For optional inclusion in 
this setup, as shown in '’PMIS 1 Program Card Setup": 

a. Tl and AT card (if such values are kn'-wn prior to program 
execution) . 

2. as a separate, single-phase job to obtain PMIS 2 data not pre~ 
viously run (where the Test, or file, has been processed by PMIS 

1, and PMIS 2 program was not included in the job setup). Optional 
inclusion in this setup, as shown in "PMIS 2 Program Card Setup": 

a. Tl and AT Card; 

b. Disk Address Card, 2110 format 


3 . as. a. separate , single-phase job to change existing PMIS .2 data 
by re-running it with new values of Tl and AT . Inclusion in this 
setup, as shown in "PMIS 2 Program Card Setup", is not optional. 

It is: 

a. mandatory that the Tl and AT card, with new, changed values 
of Tl and AT, he included in setup; 

b. mandatory that the Disk Address card, 2110 format, be included 
in the setup. 
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PHIS 2 PROGRAM CARD SETUP 


The PMIS 2 program Is executed as a separate, single-phase job by run- 
ning the following setup on the System- 370: 







" 99 ^isk ‘Address Card” - 2110 format 

1. Mandatory in PHIS 2 as single-phase program re-run with new values of T1 & AT* 

2. Optional in PMXS 2 single-phase pro gram except in above re-runs. 

3S1 . 96 ’ 3S5.S? * W j ” T1 and AT Card” 

1. Mandatory in PMIS 2 as single-phase program re-run with new values of T1 & AT. | 

2. Optional except in re-runs, as above. 
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LOG NO. . 
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@0 
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REQUESTOR ^00 p C.T 

DATE 
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.TIME. 
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PHIS 2 PROGRAM - USAGE AS SINGLE-PHASE JOB 


PHIS 2 Card Setup . Control Cards 

When the PHIS 2 program must be re-run with values of T1 and AT that 
are different from those used previously on a particular set of files 
(Tests) : 

a. the PMI3 2 program must be executed as a separate job; 

b. a Tl and AT card punched with the new values of Tl and AT must 
be included in the card setup; and* 

c. a Disk Address card punched in the 2110 format must be included 
in the card setup. 


1 . Tl and AT Card 


Card Format: 4F10.4. 

Words are: Tl vertical; AT vertical; Tl horizontal; and* 

AT horizontal, respectively. 


2. Disk Address Card 

The disk address card is used to tell the program what values of 
HIGH^FF and JljlBEND to use for this run. Unlike the Disk Address card 
used in PHIS 1, the 2110 format card used in PHIS 2 overrides the 
disk directory instead of updating it. 

Format (when used in PMIS 2 program): 2110. Words: HIGH(j5FF; 

and, J<j)BEND, respectively. 

a. The first word, HIGH^FF, is -1+ the disk address of the first file 
to be recomputed. HIGH<j>FF is the highest address which PMIS 2 

or MFMR 2 has operated on up to that time. 

b. The second word, J(j)BEND, is 149+ the disk address of the last file 
to be recomputed. JljjBEND is the highest disk address PHIS 1 

or MFMR 1 has operated on at the, time the last job was ended. 

Note: if PMIS 2 encounters MFMR data, such data are bypassed. 

Disk directory usage and the usage of HIGH(j)N and HIGH([)FF values proceed 
as if these data had been used. 
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NOTES ON CONVERSION TO ENGINEERING UNITS 
IN THE PHIS 2 PROGRAM 


HOUSEKEEPING DATA 

PHIS 2 converts the housekeeping data to engineering units by means of 
a table- look-up subroutine that performs linear interpolation between the 
closest two points. These tables are derived from tables A, B, C, D, E, and 
F for PMIS Temperature Calibration as per TM1353, Appendix K, pages 1-4. 


UNCORRECTED SKY BRIGHTNESS TEMPERATURE, T*ft 

The uncorrected sky brightness temperature is computed according to the 
formula: 

Ca “ On 

*■» " 11 + M ’ 

where: 

Ti and AT are the Y intercept and slope provided on the T^, AT card; 

(Ja is the average value of the counts when the radiometer is in the 
operate mode; 

C B is the average of the base line counts; and, 

G c is the average of the calibrate counts. 


STANDARD DEVIATION OF THE UNCORRECTED SKY BRIGHTNESS TEMPERATURE, 0 T 

The standard deviation of the uncorrected sky brightness temperature is 
computed according to the formula: 

"I ■ ' fo JqJ* c - 5s) 2 °A 2 + <?A - c<J 2 °B 2 + Pa - ^) 2 a c 2 > 
where: 

= standard deviation of data counts; 
o B = standard deviation of baseline counts; and, 

Oq = standard deviation of calibrate counts. 
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NOTES ON CONVERSION TO ENGINEERING UNITS 
IN THE PHIS 2 PROGRAM 


cont'd 

STANDARD DEVIATION OF THE UNCORRECTED SKT BRIGHTNESS TEMPERATURE, ff T 


A standard deviation of zero in any of the data is an indicator of hard- 
ware trouble. Therefore, if cr^, cr B , or a c is zero, a T is flagged by changing 
its sign. A negative value of is thus an indicator of bad data; and, since 
this sign propagates through all subsequent calculations of antenna loss and 
radozne loss, it automatically flags these calculations also. 
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SKY TEMP PROGRAM 


Before the antenna loss program (PMIS 3) can be executed, it is neces- 
sary to run program SKYTEMP with meteorological data taken on the same day 
as that of the PMIS data. 

The input data cards for SKYTEMP are of two types: (1) the header card; 

and (2) pressure, temperature, and dew point (PTD) cards. 

The Header Card, followed by the PTD cards, must fit the following for- 
mats: 


80 COLUMN PUNCH CARD FORMAT 



Card format: 3E12. 1,3X, 5A4 

The Log Id. word in the header card must be the same as the Log Id. 
word on the first data tape for that date, i.e., for the first data tape 
to be associated with the particular meteorological data. 

The sky temperature may be computed at any number of frequencies by 
specifying the start frequency, the frequency increment, and the stop 
frequency. If only one frequency is desired, ensure that stop frequency 
value is smaller than that for start frequency by leaving stop frequency 
blank (see header cards in example of SKYTEMP program setup) ; increment 
| field may also be left blank. 

; I 

2) . PTD Data Cards 

Card format: 3F10.1. 

Wherever dew point data are missing, -99.0 should be punched in 
that field of the PTD data card(s). (See PTD cards in example of SKY- 
i f TEMP program card setup.) 
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SKYTEMP PROGRAM CARD SETUP 


Any number of these data sets may be stacked with one /* card separating 
data sets. Note that two consecutive /* cards terminates the program. 

In order to run program SKYTEMP, execute the following statements on 
the System-370: 

/ kfjd of Jo h ( £ Oj) 
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APPROX RUN TIME 
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SKYTEHP PROGRAM 


Each time a sky temperature data set is computed, the program punches 
a card (the "SKYTEMP Card") which is used by PMIS 3 as input data. One SKY- 
TEMP card is punched for each sounding at each frequency requested. Normally, 
two soundings, at times spanning the time at which the radiometer measurement 
was run, are used. 

The form of the SKYTEMP card is as follows: 
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"SKYTEMP CARD" OUTPUT BY SKYTEMP FOR INPUT TO PMIS 3 PROGRAM 


Note : the "0.0" punched in card columns 51, 52, and 53, above, has no sig- 
nificance in the PMIS program (the field is still used in the MFMR program, 
however) . 
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PMIS 3 

ANTENNA LOSS FACTOR - RADOME LOSS FACTOR 


PMIS 3 is the program that computes the antenna loss factors and the ra- 
dome loss factor. 


Antenna loss is computed according to the formula: 


Lant 


W - 

T ant “ T s ky 


where T ant 


Tg + T 9 + 4 Tip 
6 


+ 273.1 °K, 


the average of the six antenna thermistor readings. 

T*jj is the uncorrected sky brightness temperature from PMIS 2 calculations. 
T s py is the sky temperature from SKYTEMP. 


Radome loss is computed according to the formula: 


LR 


T f jtR ~ Tant (1 ~ Tant) ~ Tr Lant 
TskyR T a nt “ Tr Lant 


where: 

T’BR is the uncorrected sky brightness temperature from PMIS 2 taken 
with the radome on; 

Tant is as above with radome on; and, 

Lant is the value computed above with the radome off. 

Tgfcy^ is the sky temperature from SKYTEMP at the time of the T'br 
measurement. 

Tr is the average of six radome temperatures. 


Tr 


Tn + T 12 + 4 T l3 
6 


+ 273.1 °K 
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RADOME LOSS 


PHIS 3 

STANDARD DEVIATIONS - ANTENNA LOSS, 


The standard deviation of the antenna loss is computed according to the 
formula: 


0 LA 


La cr’rj' 

S'BR ” t a 9 


where is the standard deviation of T 7 jj R . 


The standard deviation of the radome loss is computed according to the 
formula: 


J LR 


L R l A 0 T 

t'br l a - Ta (La - D - Tr 


where T R is the kinetic temperature of the radome. 
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PMIS 3 PROGRAM - ''RADOME IS OFF" CARD SETUP 



0 


0 


A PMIS 3 "Radome Is Off" job is run by executing the following program 
on the System-370: 



» First SKYTFMP 


PMIS Z Card 


JOB NAME 


REQUESTOR 


APPROX RUN TIME 


TAPE NAME & DATE CREATED 


VOLSER NO, 


TAPE LOC 


REAL ALLOCATION 


PUNCH 
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4 

COMMENTS 


SPECIAI 




CARRIAGE TAPE 
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PHIS 3 PROGRAM 


Card Setup - "Radome Is Off” 


Any number of data sets may be stacked behind the initial data control 
cards. The first data set of each job setup must use the SKYTEMP cards, but 
a /* card replaces the SKYTEMP cards in succeeding sets thus stacked. 


Data Cards 


The first data card after the // EXEC EM03C is the "Radome Is Off" 
card (card #7). 

Card i?8 is the FMIS 2 card punched by the PMIS 2 program and used 
to tell PMIS 3 where to retrieve the uncorrected sky brightness and the 
thermistor temperatures. The format of the PMIS 2 card (given in the 
PMIS 2 writeup) is repeated here for your convenience: 
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"PMIS 2 CARD" OUTPUT BY PMIS 2 POR INPUT TO PMIS 3 PROGRAM 
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PMIS 3 PROGRAM 


Data Cards, "Radome Is Off", cont'd 


3. Cards #9 and //10 are the SKYTEMP cards. One SKYTEMP card is 

punched by the SKYTEMP program for each sounding at each frequency re- 
quested. Normally, two soundings, at times spanning the time at which 
the radiometer measurement was. run, are used. If only one sounding Is 
used, the second SKYTEMP card is left blank. Note that two cards are 
required, i.e., a blank card must be used to represent the second SKY- 
TEMP card. 

The sounding must be taken on the same day as the radiometer 
measurement, and the LOG ID on the SKYTEMP card must be repunched to 
match the LOG ID on the firs-. PMIS 2 card. Failure to observe these 
rules will produce unpredictable results. 

The format of the SKYTEMP card (given in the SKYTEMP writeup) is 
repeated here for your convenience: 


4J /*> 


3 

& § 
cn N 


a) 

4-1 

Q 


OJ 

n 

E-» 


TJ 

1-1 

bO 

o 

-3 


!>. 

u 

g 

o* 

41 

£ 



Ofl 00 0 0 0 0 0 DOfl 0 0 00 0 0 0 0 0 0 0 DO 0 0 0 0 0 0 0 0 0 0 0 0 : 0 0 0 0 ! 0 I 0 0 0 0 0 I* 0 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
i j 3 « s i 7 a 9 rn n a a it s ic 11 is n nn n n jus zs nnisjo 31 u nu jsk jnias «<i« *h(<s*'«io«j 53 Si 5 ? 5 :w 5 SK 5 is. -- si dmgieg 575353 jo n ran it ts ts n rattsa 
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555555555 555555 55555555555555555555555555 555555555555555555555555555555555555 
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"SKYTEMP CARD" OUTPUT BY SKYTEMP FOR INPUT TO PHIS 3 PROGRAM 


Each time PMIS 3 is executed for a "Radome Is Off" job, the program pun- 
ches a "Radome Is On" card for use later when a PMIS 3 "Radome Is On" job is 
run. 
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PMIS 3 PROGRAM - "RADOME IS ON" CARD SETUP 


Any number of these data sets may be stacked with one /* card separating 
data sets. 

To execute PMIS 3 when the "Radome Is On", run the following program: 


' BOJ Cent* Job') 
EOF (e*et ef 



Hi Sftfc C-ar-d. ir] jofatie o-f- SaG&ncL •S/fY'f MfrfP card 
6*70. ^,curud * *7% *s be-in $ ti-se-d h^f t ’<z) 

os #. o F/r ^ t s#yt£mp card 

— — «i — i — - — 1 * *' — 


BE C. £ 

F’-'F?: FMO'ir 

// fy 70r i , ro- nni , ] . p, nF.nors, ftp ? 00 

// n: St. F-IAtfF* ’niim 

// Pf>M?r fjniiwr rPZOOJ HW ??5 

/./ iof rnnPFR-^i oo-p» 1 4305 : pwi-s 


ii . u _.. . uu. tiu.*.. 


LOG NO. 

REQUESTOR. 
OATF 


(APE DR 


JOB NAME "Pf^SXS 3 

W. Coo^g^ TEL. 383 

TIME APPROX RUN TIME 


TAPE NAME & DATE CREATED 


®D FUND t*/3GS W n 'PMt.S 3 
CC #./ Of? RPniN 


3 


PUNCH. 


2701 SEL. 


££* BEGIN 

f wn 


TAPE LOC. |) DISK VOLSER NO. 


33S* WSOOl 




REAL ALLOCATION 


FORMS 1 2 4 

SPFCtAI 

COMMENTS 


CARRIAGE TAPE 




1 


950 Rev. 8-74 


£2*^ PASS «, 
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PMIS 3 PROGRAM 


Card Setup, "Radome Is On 11 

The same card setup is used as for the "Radome Is Off" job, except that 
the "Radome Is On" card will now be the f: rst data card in each data set. 

Any number of data sets may be stacked behind the initial data control cards. 
The first data set of each job setup must use the SKYTEMP cards, but a /* 
card replaces the SKYTEMP cards in succeeding date? sets. Radome on and off 
sets may be intermixed in the same job deck. 

The LOG ID of the first PMIS 2 card is checked against the LOG ID on the 
SKYTEMP cards. 


"Radome Is On" Data Card 


The "Radome Is On" card contains the antenna loss values needed when 
computing the radome loss. 

' Rfflg mE’T & TI Fr M-DS5V* ti.i M 

• * ' “ i,. , . - .V- ' t v| 1 Of-’? 

f. ■ .%• -v, J.. ■ 


0 0 8 0 0 0 0 0 0 0 0 04J 0 O’JO ff(f ,0 0 

...oil m ; m s id ir a lin istei) v 

» t A 4 m a ■ * * * 4 '_+] -* '«_4 


'I . 


0 0 0 0000 0 0 D 00; 

I|3Q,3I353MJ535 3) 381S4N1 « 43 « 4Hfr 

: ~ y ~ i t 1 1 1 1 U 11 lli 


This is an average of the ALOSSV 
iM values taken from several "Radome \i22 2 222 2 2 2 ; | 


|- On" cards. 


*■ . « *J , 


^ ^ * • " » i 

•SS'OS'S^St 1 




rc-J3 3 3 39^ 3 3 3 3 


■4 


0 0‘ 0 0 0 0 QOODQOOOa'0'l 

Si ibj h mss sf si sno sj si e 5 M ties efti 

|1 1 1 1=1 1 } 1 11 1 1 1 1,1 Tf 1,1 

fl 2 2 2 2 2 2 ’2 ? 2 2 2 2^ 

3 3 ai'hll’lt? 3 3 U 3 
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This is an average of AL05SH 
the same cards used to obtain the average 

ALOSSV. | l~ | 
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i . . 'V , * - . ' . - i .* 1 * rt ■■ ■ - " ' i ' 1-1 ■ * ■ iiii, ~' — r * 1 H 


• - * - . y - 

values taken from * 4 4.4 

! f " 

,5«55 



'5 S 5$ 


1 

’I'll 1'1 

.iSti ! 
uiisfP 

J|gg 


= 'Ml 1111^2 liillTlf'tlii This the disk address taken from one of the same cards 

used to obtain the average values of ALOSSV and ALOSSH. 
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1 tiw >1 
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tl 1.4 1 B J 1 l U n ti [tn 15 ts n BI3J371 IlaUKlfrajlrtW 513! 3334 M3#31 JUJ tOtt «IW« 45t6 4H34JH 5t S9i3 ttt? PUCS SW 8863 R n.IUi H lire n fatiS * s/ 
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"RADOME IS ON" CARD - OUTPUT BY PMIS 3 "RADOME IS OFP" JOB 


Values punched in the "Radome Is On" card are: antenna loss vertical, 

antenna loss horizontal, and, disk address, respectively. 

The disk address on the "Radome Is On" card is used only to check the 
beam position, and it is therefore possible to use the average of several an- 
tenna losses on this card, giving the disk address of one. 
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MFMR PROGRAMS - REFERENCE INDEX 

Page No.’s 

GENERAL M.l - M.5 


A. To load and execute System-7 PCM Program M.l 

B. Instructions of A., above, for newcomer to System-7 ..... M.la 

C. Telemetry Playback Console - sketch and power settings. ... M.2 

D. Comments on System-7 Program. M. 3 

E. Example: printout from control typewriter, MFMR data run . . M.4 

F. Notes on control typewriter printout 

1. Explanation of terms M.5 

2. "Overflow". M.4, 5 

3. "Valid-Data Flag-Drop". M.4, 5 

MFMR 1 PROGRAM Ml.l - Ml. 3 


A. Card setup for running MFMR 1 on computer Ml.l 

B. Discussion of MFMR 1 card setup 

1. Phase 1 (MFMR 1) Ml. 2 

2. Phase 2 (MFMR 2) Ml. 2 

3. Phase 3 (Job-Step Disk-Dump) Ml. 2 

C. Optional Control Cards - MFMR 1 (1st phase) 

1. UPSI card Ml. 2 

2. Disk address card (as used in MFMR 1) .... Ml. 3 

D. Optional Control Cards - MFMR 2 (2nd phase) 

1. // EXEC EM03E card, /* card Ml. 3 

MFMR 2 PROGRAM M2.1 - M2. 4 


A. MFMR 2 Program as 2nd phase of a MFMR 1 job setup 

1. "MFMR 2 Card" output by MFMR 2, input to MFMR 3 M2.1 

2. Job-Step Disk-Dump M2. 2 

B. MFMR 2 Program as single-phase, separate job 

1. MFMR 2: separate job or 2nd phase treatment M2. 3 

2. Disk address card: when used, when omitted ....... M2. 3 

3. Disk address card (as used in MFMR 2): purpose, format. . M2. 3 

C. Card setup for running MFMR 2 (separate job) on computer. . . M2. 4 

SKYTMFMR PROGRAM MS.l - MS. 3 


A. Input data cards for SKYTMFMR program 

1. Header card, description and format MS.l 

2. Pressure, temperature, dew point (PTD) cards - format . . MS.l 

B. Card setup for running SKYTMFMR on computer MS. 2 

C. "SKYTEM? Cards" output by SKYTMFMR, input to MFMR 3 MS. 3 
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MFF1R PROGRAMS - REFERENCE INDEX, confd 


Page No . 1 s 

MFMR 3 PROGRAM M3.1 - M3. 9 


A. Uncorrected sky brightness temperature, T’b M3.1 

B. Standard deviation of uncorrected sky brightness 

temperature, Otb . M3.1 

C. Antenna loss . M3. 2 

D. Standard deviation of antenna loss. M3. 2 

E. Radome loss . M3. 2 

F. Standard deviation of radome loss M3. 2 

G. Card setup for running MFMR 3 "Radome Is Off" job 

on computer M3, 3 

H. Discussion: card setup for "Radome Is Off" job deck 

1. Stacking data decks M3. 4 

I. Control cards for "Radome Is Off" job deck 

1. "Radome Is Off" card M3. 4 

2. "MFMR 2 Cards": OPR, CAL, and BASE M3. 4 

3. T1 card M3. 5 

4. AT card M3. 6 

5. "SKYTEMF Cards" M3. 7 

J. Card setup for running MFMR 3 "Radome Is On" job 

on computer M3. 8 

K. Discussion: card setup for "Radome Is On" job deck 

1. Stacking data decks . M3. 9 

L. Control Cards: "Radome Is On" card M3. 9 
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PROGRAM PROCEDURE TO REDUCE MFMR DATA TAPES 

PRIOR TO LOADING THE SYSTEM-7 PCM PROGRAM! 

PERFORM CHECKOUT OF ALL PCM HARDWARE AND TAPECS) , 

TO LOAD AND EXECUTE THE SYSTEM-7 PROGRAM; PROCEED AS FOLLOWS: 

(1) RUN U-ZERO TAPE (LOAD UNIT ADDRESS = X'0000') 

(2) LOAD IPL PROGRAM FROM DISK (LOAD UNIT ADDRESS = X'0002') 

*** I POO A ENTER CONTROL STATEMENT 

(3) CYCLE HOST ATTACH SWITCH TO (ENABLE AND IPL) AND 

THEN TO (ENABLE).i 

IT MUST REMAIN ON (ENABLE) FOR REMAINDER OF PROGRAM EXECUTION. 

tt) SET DISK PROGRAM "DATA SET NAME" BY FOLLOWING R(EFER) 
STATEMENT: 

R JOBLIB;F2;;PCMUSER 

(5) L(OAD) PCM PROGRAM INTO CORE BY FOLLOWING L(OAD) STATEMENT: 

L PCM01B 


THE SYSTEM-7 IS NOW PREPARED TO ACCEPT LINKUP WITH THE SYSTEM-370 
AND WILL REMAIN IN A WAIT STATE UNTIL MFMR 1 PROGRAM IS EXECUTED 
BY THE SYSTEM-370. 
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PROGRAM PROCEDURE TO REDUCE MFMR DATA TAPES 


•'■ cites : "MFMR" is the name of the program and of the radiometer being tested; procedures 

are given in large type, while explanatory comments are given in smaller type enclosed 
by hand-drawn brackets, thus: [For benefit of one who is unfamiliar with the System-7] . 

PRIOR TO LOADING THE SYSTEM-7 PCM PROGRAM: 

PERFORM CHECKOUT OF ALL PCM HARDWARE AND TAPE(S). 


TO LOAD AND EXECUTE THE SYSTEM-7 PROGRAM, PROCEED AS FOLLOWS: 

(1) RUN U-ZERO TAPE (LOAD UNIT ADDRESS = X'0000') 

|u-ZERO, the tape to IPL the computer, is a small, blue punch tape kept 
right by the computer; put Address switches on OOOO (Addresses are con- 
trol switches on front of computer; there are 4 Address switches, thus 
each Address switch is set on 0, as indicated by the four zeros in 
step (1), above], 

(2) LOAD IPL PROGRAM FROM DISK (LOAD UNIT ADDRESS = X'0002') 

[iPL denotes "Initial Program Load", Leave the first 3 Address switches 
on zero, move the right hand Address switch to setting of Q 

*** I POO A ENTER CONTROL STATEMENT 

jthe above is a statement that the computer sends back to you after you 
have done step (2)]] 

(3) CYCLE HOST ATTACH SWITCH TO (ENABLE AND IPL) AND THEN TO (ENABLE) 
WHERE IT MUST REMAIN FOR THE REMAINDER OF PROGRAM EXECUTION, 

[the Host attach switch is another contrc switch on the front of the com- 
puter; switch it to "Off" (that’s down), then up to "Enable and IPL" and 
back to center position, which is "Enable"; it must stay on "Enable" all 
the time you’re transferring data from Sys.-7 to Sys.-370 - this is done 
in order to establish a linkage with the System-370T 

m SET DISK PROGRAM "DATA SET NAME" BY FOLLOWING R(EFER) STATEMENT: 

[step (4) tells you to type an "R", that is, "Reference Statement", on the 
control typewriter. The statement to be typed, is as follows, below (be 
sure to leave a space between the "R" and the 

R J0BLIB,F2,,PCMUSER 

(5) L(OAD) PCM PROGRAM INTO CORE BY FOLLOWING L(OAD) STATEMENT: 

[after typing the statement of step (4)’s instruction, execute step (5) by 
typing the following statement, below. Leave a space between "L" and "P’Q 

L PCM01B 

THE SYSTEM-7 IS NOW PREPARED TO ACCEPT LINKUP WITH THE SYSTEM-370 AND 
WILL REMAIN IN A WAIT STATE UNTIL MFMR 1 PROGRAM IS EXECUTED BY THE 
SYSTEM-370. 
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COWOTS 0^1 SYSTEM-7 PROGRAM page i 


oTMT 


SOURCE STATEMENT DOS ASM/7 

(360A— TX— Oil) V1M2 

12/10/74 

2 

ft 

PCN010 PROGRAM SPECIFICATIONS: 




00000020 

•3 

ft 

PURPOSE: INPUT { MFMR ) OR (PM IS) 

DATA INTO A 649 WORD BUFFERS 

00000040 

4 

ft 

(A GR B), DETERMINE TYPE { MF MR ) OR (PMIS) 

f CHECK FOR 

00000050 

5 

£ 

PROPER SEQUENCE OF (TIME-WORD- 

3-FLAG) AND 

(FID) FORMS f 

00000060 

6 

* 

TRANSFER TOGGLED BUFFER 

AT (TIME-WORD-3) 

IF (VALID DATA) 

00000070 

7 

fit 

FLAG IS PRESFNT, AND SIGNAL ( ENO-OF-FILEI 

WHEN 

OOOOOOGO 

e 

ft 

(VALIO-DATA-FLAG) DROPS, 




00000090 

9 

ft 





000O0L00 

10 

•ft 





00000110 

11 

ft 

BUFFER FORMAT; 




ooooouo 

12 

ft 





00000 L 30 

12 

ft 

HEADER (8-WORDS) DISP 

(0000) 

CODE WORD 

(0000) DATA 

00000140 

14 

ft 




t FFFF) EOF 

00000 150 

15 

ft 

o © 
**1 

t=! 

o| 



( EEEE) EOJ 

00000160 

16 

ft 


BLANK 


00000170 

17 

* 

(0002) 

LOG ID 2-CHARS (EBCDIC) 

00000180 

18 

ft 

8| 

(0003) 

LOGID 2-CHARS (EBCDIC) 

00000190 

19 

ft 

(0004) 

TIME-WORDI 


00000200 

20 

ft 

(0005) 

TIMF-WQRD2 


00000210 

21 

ft 

^ § 

(0006) 

TIME-WORD? 

(QUEUE XFR ) 

00000220 

22 

ft 

S H 

(0007) 

BLANK 


00000230 

23 

*ft 

i — > 
*3 CO 




00000240 

24 

ft 

FRAME NUMBER tl) 

(0008) 

SYNC1 


00000250 

25 

ft 


(0009) 

SYNC 2 


00000260 

26 

ft 

DAT Al 1 (30 WORDS) 

(0010 ) 

DATA 

• 

00000270 

27 

ft 


TH DI J 



00000280 

28 

ft 


(0039) 

DATA 


00000290 

29 

ft 





00000300 

30 

ft 

F ID I 

(0040) 

DID1 (FID 

CHK j SET TYPE) 

00000310 

31 

ft 





00000320 

22 

ft 

DATA12 ( 127 WOROS) 

(0041) 

DATA 


00000330 

33 

ft 


THRU 



00000340 

34 

ft 


(0167) 

DATA 


0G000350 

25 

ft 





00000360 

36 

ft 

FRAME NUMBER (2) 

(0160) 

THRU (0169) 

SYNC (TIME CHK J 

00000370 


GO 

fD 


s 

U> 


EXAMPLE OF PRINTOUT FROM CONTROL TYPEWRITER - MFMR DATA RUN 


The following System-7 messages will be executed by the System-7 when the Sys- 
tem-370/System-7 linkup occurs, as indicated by the statement ”S370 is ready 
to accept data”. Responses typed by the operator on the control typewriter 
are indicated by enclosure in a box. 


*** IFOO A ENTER CONTROL STATEMENT 
j R JOBLIB.F2, .PCMUSER I 
I L PCM01B " 

S370 IS READY TO ACCEPT DATA 

ENTER OPTIONS BY TYPEWRITER (REQUEST) : 

REQ (LOG) TO ENTER LOG-ID WORD THEN (OSC) OR (RUN). 

REQ (OSC) TO RUN SIMULATED PCM DATA. 

REQ (RUN) TO RUN PCM DATA. 

REQ (EOF) MARK PILE AND IGNORE DATA. THIS IS NORMALLY 
AUTOMATIC BY VALID-DATA-FLAG DROP. 

REQ (EOJ) TO SIGNAL TERMINATE JOB. 

NOTE — NEW FILES MAY BE OPENED BY (LOG) /(OSC) OR (RUN) REQUESTS 
AFTER AN (EOF) REQUEST. 


LOG 


ENTER 4-DIGIT (HEX) LOGBOOK(ID) WORD 


I MQ19 

REQUEST (OSC) OR (RUN) TO PREP FOR PCM DATA 


OR 1; | RUN 


ENTER PCM DATA 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 
EOF PROCEDURES 


EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 

EXECUTED 


overflow occurred 


fsee following page, ’’Notes on 
J Control Typewriter Printout - 
MFMR Data Run” 


l 


OR 1:1 RUN 


ENTER PCM DATA 

EOF PROCEDURES EXECUTED 

EOF PROCEDURES EXECUTED 


job was started again at selected place on tape 


ORl: [ EOJ~| — — z — — — • end of last tape; and of job was requested 

S7/S370 LINK TERMINATED 
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NOTES ON CONTROL TYPEWRITER PRINTOUT - MFMR DATA RUN 


Explanation of Terms 

LOG ID = name of the tape (4-digit word) . The first digit must be 
either "M" or "P" in order for the computer to recognize 
data as being either MFMR or PMIS data. 

OSC - running oscillator to de-bug tape. 

PCM = pulse code modulation. 

VDC = Test Id. number (name of file). The Test Id. is the 
"valid data" code, or word. 

EOF = end of file. 

EOJ - end of job. 


Overflow 


When overflow occurs, as it did during the processing of MFMR tape 
M019 (indicated by arrow on "Example of Printout from Control Typewriter - 
MFMR Data Run" on the previous page), determination of where to re-start 
the job is made by: counting EOF f s; and, by checking the time on the time 
code reader and the time code on the log. The tape is then backed up to 
the selected point, and the job is started again. 

Although whole files are not lost prior to overflow (all files are lost 
afterward), it is possible to lose 25 data cycles, spread throughout the 
files, out of the job. There are 640 words per data cycle. 

When data overflows: a belJ rings; EOF messages cease to print out; 

and, the control typewriter skips lines where the messages would have been. 


Valid-Data-Flag Drop 

The twenty-second word of each PCM data cycle is the "valid data word" 
(VDC). The high-order bit of this word is the "Valid-Data Flag". The 
valld-data-flag is one (1) during each file and is zero (0) between files. 

The System-7 program monitors the valid-data-flag and executes an end- 
of-file procedure each time the valid-data-f larv ’rops, i.e., goes from 1 to 0. 
The MFMR 1 and PMIS 1 programs also monitor the ✓alue of the entire valid- 
data-word; in case the System-7 fails to execute an end-of-file procedure 
any time the VDC changes, the MFMR 1 and PMIS 1 programs will do so. 
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The MFMR 1 program is executed by running the following setup on the 


(MFMFi 1 ) 


APPROX RUN TIME 


TAPE NAME & DATE CREATED 


TAPE LOC. 


VOLSER NO, 


TAPE DR 


DISK 


















MFMR 1 PROGRAM CARD SETUP & PROGRAM OPTIONS 


HFMR 1 Program Card Setup 

1. Phase 1 

The example of a "MFMR 1 Program Card Setup" given on the previous 
page actually includes three separate programs, MFMR 1, MFMR 2, and Job" 
Step Disk-Dump, all executed by the MFMR 1 program setup. When the three 
programs are thus used, MFMR 1, MFMR 2, and Job-Step Disk-Dump are the 
first, second, and third phases, respectively, of a single job. 

By removing the MFMR 2 cards numbered 12 and 13, MFMR 1 can be run 
as a separate job. 

2. Phase 2 

MFMR 2 may also be run as a separate job. 

3. Job-Step Disk-Dump 

When MFMR 2 (or MFMR 1, if MFMR 2 was not included the job setup) 
has processed all of the data on the disk, control is turned over to Job- 
Step Disk-Dump (cards #14 and #15). Job-Step Disk-Bump prints all of 
the data cycles that were dumped into the disk by MFMR 1. 

This job step also computes the amount of run time for the entire 
program and prints this time on the SYSLOG and on the printer, thus must 
be included in all foreground MFMR and PMXS programs. Disk Dump dumps 
the disk only if it is enabled by following the // EXEC EM03H card (#14) 
with the dump card (#15). If no dump is desired, follow the // EXEC 
EMG3H with the /* card (#16). 


Optional Control Cards - MFMR 1 (First Phase) 

1, Card No. 8, the UPSI card, is optional. It contains a string of 

eight binary digits. 

The binary value of the seven low-order bits of this string is the 
number of data cycles which will be dumped onto the disk each time a new 
Test Id. is established, l.e., the number of data cycles dumped per file. 
Note: the last data cycle of each file (Test) is not under the control 

of the // UPSI card and Is always placed (dumped) on the disk. 

The high-order bit is a flag to determine whether or not the diag- 
nostic Is to he printed from MFMR 1 (Phase 1 of the example, "MFMR 1 
Program Card Setup"). If the high-order UPSI bit is 1, the full print- 
out is obtained; if the high-order bit is zero, or if there is no UPSI 
card. Phase 1 printout is suppressed. 
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MR 1 PROGRAM CARD SETUP & 


PROGRAM OPTIONS 


Optional Control Cards - MFMR 1 (First Phase), contM 

2. Card 10, "Disk Address" card, must be used with great caution. 

Its misuse can destroy all data on the disk, since it updates, i.e., 
changes the disk directory. On the card shown in the MFMR 1 program 
setup, each word has the value n 99", which is the starting address. 

This initialises the di3k and causes it to start loading data at the 
beginning of the disk file. 

Format of Disk Address Card (when used in MFMR 1) : 3110. 

a. First word, HIGHER: the highest disk address that MFMR 1 or 

PMIS 1 has operated on up to a given instant; 

b. Second word, HIGH(j)FF: the highest disk address that MFMR 2 or 

PMIS 2 has operated on up to the same given instant; 

c. Third word, J({)BEND: the highest disk address that the last 

MFMR 1 or PMIS 1 job has operated on up to the time it was 
terminated. 

Note: the disk control is set up such that both MFMR 1 and MFMR 2 

can operate simultaneously in different partitions of the computer. 

And, if MFMR 1 is in process of executing while MFMR 2 is also in the 
process of executing, then MFMR 2 has to know exactly where MFMR 1 is. 
Three variables are transmitted through the disk so that these two 
programs can communicate with one another; these variables are HIGH(j>N, 
HIGH(|>FF, and J(j>BEND. 


Optional Control Cards - MFMR 2 (Second Phase) 


3. Cards numbered 12 and 13, the // EXEC EM03E and the /* cards, 

respectively, are optional in the MFMR 1 program. Inclusion of these 
cards causes the MFMR 2 program to be executed In the MFMR 1 program 
job (see example of MFMR 1 program card setup). 

When the MFMR 2 program has not been included in a MFMR 1 job and 
must be run at a later time, MFMR 2 is executed as a separate job (ref. 
"MFMR 2 Program - Usage as a Single-Phase Job"), 

If it becomes necessary to re-run MFMR 2 data that were run previ- 
ously, the MFMR 2 program must be executed as a separate job (ref. 

"MFMR 2 Program - Usage as a Single-Phase Job") . 
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MR 2 PROGRAM 


When the EOJ instruction is issued at the System-7 terminal, the link- 
age between the System-370 and the System-7 is terminated, and program MFMR 1 
(the first phase of the job) ends and turns control over to the second phase 
(MFMR 2 program). 

MFMR 2 retrieves the raw data placed on disk DRZQ01 by MFMR X and com- 
putes the engineering units; it then places these results on the disk, prints 
them out, and punches a card for each file (Test). These cards are referred 
to as "MFMR 2 Cards" (see example) . 

The MFMR 2 card will be used in the MFMR 3 program (Ejee example of 
MFMR 3 program card setup) . 
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"MFMR 2 CARD" OUTPUT BY MFMR 2 PROGRAM 


99! 


Mote: the Log Id., above, is the name of the tape. The first digit must be 
either "M" or "P" in order for the computer to recognize data as being either 
MFMR or PMIS data. The Valid Data Word (or Valid Data Code, VDC) is the Test 
Id. number (name of file). 
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MFMR 2 PROGRAM 


VIhen MFMR 2 has processed all of the data on the disk, control Is turned 
over to Job-Step Disk-Dump. This is the last job step of the three-phase job 
illustrated by the MMR 1 Program Card Setup. Job-Step Disk-Dump prints all 
of the data cycles that were dumped into the disk by MFMR 1. Note: the last 
data cycle of each file (Test) is not under the control of the / / UPSI card 
and is always placed (dumped) on the disk. 

Job— Step Disk— Dump also computes the amount of run time for the entire 
program and prints this time on the SYSLOG and on the printer „ thus must be 
included in all foreground MFMR and PMIS programs. Disk Dump dumps the disk 
only if it is enabled by following the // FiXEC EM03H card with the dump card. 
If no dump is desired, follow the // EXEC EM03H with the /* card. 
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MFMR 2 PROGRAM - USAGE AS SINGLE-PHASE JOB 


MFMR 2 Card Setup. Single-Phase Job . 

The MFMR 2 program may or may not be run as the second phase of the three- 
phase MFMR 1 job as shown in the MFMR 1 program card setup; its inclusion is 
optional. However, if the MFMR 2 data are to be obtained at a later time, or 
if it should become necessary to re-run MFMR 2 data, the MFMR 2 program must 
be executed as a separate, single-phase job. 


1, Run MFMR 2 as a separate, single-phase job to obtain MFMR 2 data 
not previously run (where the Test, or file, has been processed by 

MFMR 1, and MFMR 2 program was not included in the job setup) as follows: 

a. with card #7, "Disk Address" card omitted , run the job shown in 
the "MFMR 2 Program Card Setup". 

2. Run MFMR 2 as a separate, single-phase job to change existing 
MFMR 2 data by re-running it as follows: 

a. run the job shown in the "MFMR 2 Program Card Setup"; and, 

b. include card #7, "Disk Address" card, 2110 format. 


Disk Address Card 

The disk address card is used to tell the program what values of HIGH<j)FF 
and J<j)LEND to use for this run. Unlike the Disk Address card used in MFMR 1 , 
the 2110 format card used in MFMR 2 overrides the disk directory instead of 
updating it. 

Format (when used in MFMR 2 program): 2110. Words: HIGHljJFF, and, 

J<|iBEND, respectively. 

a. The first word, HIGHOFF, is -1+ the disk address of the first file to 
be recomputed. HIGHOFF is the highest address which MFMR 2 or PMIS 2 
has operated on up to that time. 

b. The second word, JtjiBEND, is 149+ the disk address of the last file to 
be recomputed, J(j>BEND is the highest disk address MFMR 1 or PMiS 1 
has operated on at the time the last job was ended. 

Note: if MFMR 2 encounters PMIS data, such data are bypassed. Disk di- 

rectory usage and the usage of HIGH(j>N and HIGHOFF values proceed as If these 
data had been used. 
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MFMR 2 PROGRAM CARD SETUP 
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The MFMR 2 program is executed as a separate, single-phase job by run- 
ning the following setup on the System-370: 



REQUESTOR TEL >£ & 3 CC #11 1£ BEGIN 

DATE TIME APPROX RUN TIME END _ 


TAPE DR. 

TAPE NAME & DATE CREATED 


TAPE LOG. 

DISK 

VOLSER NO. 

280 





2>f?ZOO) 

281 






282 






283 







PUNCH 270) SEL. , REAL ALLOCATION 


FORMS 1 2 

4 

COMMENTS 


SPECIAI 

CARRIAGE TAPE _ 
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MFMR 2 PROGRAM - SETUP FOR SEPARATE, SIHGLE-PHASE JOB 
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SKYTHFMR PROGRAM 


Before the loss program (MFMR 3) can be executed, it is necessary to 
run program SKYTMFMR with meteorological data taken on the same day as that 
of the MFMR data. 

The input data cards for SKITMFMR are of two types: (1) the header 

card; and (2) pressure, temperature, and dew point (PTD) cards. 

The header card, followed by the PTD cards, must fit the following for- 
mats: 


80 COLUMN PUNCH CARD FORMAT 



CARD FORMATS - HEADER CARD, PTD CARDS 

1) . SKTTMFMR Header Card 

Card format: 39X, 5A4 

The Log Id. word in the header card must be the same as the Log Id. 
word on the first data tape for that date, i.e., for the first data tape 
to be associated with the particular meteorological data. 

2) . PTD Data Cards 

Card format: 3F10.1. 

Wherever dew point data are missing, -99.0 should be punched in 
that field of the PTD data card(s). (See PTD cards in example of 
SKTTMFMR program card setup.) 
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SKYTMFMR PROGRAM CARP SETUP 


Any number of these data sets may be stacked with one /* card separating 
data sets. Note that two consecutive /* cards terminates the program. 

In order to run program SKYTMFMR, execute the following statements on 
the System-370: 

' Ft* d. & £ t> ( Jr O J ) 

* o f Fjtza tr> r/i/s (£VF) 



PUNCH 


FORMS 1 2 4 

SPECIAL 

CARRIAGE TAPE 


COMMENTS 


950 Rev. 8-74 
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SKYTMFF1R PROGRAM 


Each time a sky temperature data set is computed, the program punches 
four cards (the "SKYTEMP Cards") which are used by MFMR 3 as input data. 

One SKYTEMP card is punched for each sounding at each frequency requested. 

Following is an example of the card format and the four SKYTEMP cards; 
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"SKYTEMP CARDS" OUTPUT BY SKYTMFMR FOR INPUT TO MFMR 3 PROGRAM 
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MFMR 3 

UNCORRECTED SKY BRIGHTNESS TEMPERATURE 
ANTENNA LOSS FACTOR - RADOME LOSS FACTOR 


MFMR 3 is the program that computes the uncorrected sky brightness tem- 
perature, the antenna loss factors, and the radome loss factor. 


UNCORRECTED SKY BRIGHTNESS TEMPERATURE, T'g 


The uncorrected sky brightness temperature is computed according to the 
formula: 

T'b - Ti + AT / , 

C C “ Cg 

where: 

Tj and AT are the Y intercept and slope provided on the Tj and AT cards; 

T 3 A is the average value of the counts when the radiometer is in the 
operate mode; 

C33 is the average of the base line counts; and, 

<Tc is the average of the calibrate counts. 

STANDARD DEVIATION GF THE UNCORRECTED SKY BRIGHTNESS TEMPERATURE, cr TE 


The standard deviation of the uncorrected sky brightness temperature is 
computed according to the formula: 

| | j ... — ■ ■ — ■ . . - - ■ — 

"'TB = 7 ” ZTz G - C’bJ a A 2 + C^A - Cc) 2 (J3 2 + (Ga - C B ) 2 0 C 2 , 

(G C - C B ) 2 

where: 

cfa = standard deviation of data counts; 

(Tg = standard deviation of baseline counts; and, 

<Jc = standard deviation of calibrate counts. 


A standard deviation of 2ero in any of the data is an indicator of hard- 
ware trouble. Therefore, if or a C is zero, a-p is flagged by changing 

its sign. A negative value of is thus an. indicator of bad data; and, since 
this sign propagates through all subsequent calculations of antenna loss and 
radome loss, it automatically flags these calculations also. 
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MFMR 3 

ANTENNA LOSS FACTOR - RADOHE LOSS FACTOR 


AMTEKMA LOSS 

Antenna loss is computed according to the formulas 

_ Tg — ^A 

A " T’b Lw - *w (Ly - 1) - Ta • 

where: 

Tg is the sky temperature in °K; 

Ta is the kinetic antenna temperature in °K; 

T'b is the uncorrected sky brightness temperature in °K; 
Ly is the wave guide loss; and, 

Ty is the wave guide temperature in °K. 


STANDARD DEVIATION OF ANTENNA LOSS 


The standard deviation of the antenna loss Is computed according to 
the formula: 


L a L w o tb 

0 LA t’ b Ly + Ty (1 - Ly) - T A 


RADOME LOSS 


The radome loss is computed according to the formula: 


Lr 


?s - t R 

t'b l a l w - t a (l a - i> - Tr - t w (l w - l) l a 


where T^ is the temperature of the radome °K (average of 5 radome thermistors). 


STANDARD DEVIATION OF RADOME LOSS 


The standard deviation of the radome loss is computed according to the 
formula: 


°LR 


~ l R l A q LA 

T’b XA l W " T A (LA - 1) - tr 
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MFMR 3 PROGRAM - "RADOME IS OFF" CARD SETUP 


A MFMR 3 "Radome Is Off" job is run by executing the following program 
on the System-370: 


kCJ (Cfjti ef Jrb) 


as 





First data set must be complete, but succeeding sets may be truncated in any of 
3 places: 1) after MFMR 2-OPR card, as in 2nd data set; 2) after MFMR 2-BASE 
card; or, after AT card. 


£hll)utt 
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BEGIN 
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log no 
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RFOUESTOR 


W/* 7 . Coot 


APPROX RON TIME 



PUNCH 
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2701 SEL FIEAL ALLOCATION 


FORMS 1 2 4 COMMENTS 

SPECIAL 

CARRIAGE TAPE 


ORIGINAL PAGE IS 
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MFMR 3 PROGRAM 


Card Setup - "Radome Is Off" 

Any number of ^ata sets may be stacked behind the initial data control 
cards. The first data set of each job setup must use: the OPR, CAL, and 

BASE cards; the T1 and AT cards; and, the SKYTEMP cards. A /* card may be 
used to truncate succeeding data sets in any one of three places: after the 
OPR card; after the BASE card; or, after the AT card. 


Data Cards 


1. 


2 . 


The first data card aft«r the // EXEC EM03F is the ’’Radome Is Off” 
card (card Hi). 

Cards it 8, //9, //10 are the MFMR 2 OPR, CAL, and BASE cards punched 
by the MFMR 2 program and used to tell MFMR 3 where to retrieve the engi- 
neering data. The format (given in less detail in the MFMR 2 writeup) 
and an example of each of the MFMR 2 cards is shown below: 
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"MFMR 2 CARDS" OUTPUT BY MFMR 2 FOR INPUT TO MFMR 3 PROGRAM 
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MFMR 3 PROGRAM 


Data Cards, "Radome Is Off", cont'd 

3. Cards #11 and #12 are the Tl and AT cards required for the uncor 

rected sky brightness calculations. 


Formats 

for 

these cards 

are as follows: 
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Tl CARD FOR INPUT TO MFMR 3 PROGRAM 
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MFMR 3 PROGRAM 


I) 


Data Cards, "Radome Is Off", cont'd 
3. continued 

Format for the AT card is as follows: 
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AT CARD FOR INPUT TO MFMR 3 PROGRAM 
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MR 3 PROGRAM 


Data Cards, "Radome Is Off 1 *, cont*d 

4, Cards #13, #14, #15, and #16 are the SKYTEMP cards. One SKYTEMP 

card is punched by the SKYTMFMR program for each sounding at each of the 
four MFMR frequencies. Note that four cards are required: one card for 

each frequency in the correct order, i.e., 18, 37, 22, and 1.4 GHz, re- 
spectively. 

The sounding must be taken on the same day as the radiometer meas- 
urement, and the Log Id. on the SKYTEMP cards must be repunched to match 
the Log Id. on the first MFMR 2 card. Failure to observe these rules 
will produce unpredictable results. 

The format of the SKYTEMP cards (given in the SKYTEMP writeup) is 
repeated here for your convenience: 
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"SKYTEMP CARDS" OUTPUT BY SKYTMFMR FOR INPUT TO MFMR 3 PROGRAM 
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Each time MFMR 3 is executed for a "Radome Is Off" job, the program pun- 
ches a "Radome Is On" card for use later when a MFMR 3 "Radome Is On" job is 
run. 
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jvjpMR 3 PROGRAM - "RADOME IS ON" CARD SETUP 
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A MFMR 3 "Radome Is On" job is run by executing the following program 
on the System-370: 
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MFMR 3 PROGRAM 


Card Setup, "Radome Is On" 

The same card setup is used as for the "Radome Is Off" job, except that 
the "Radome Is On" card will now be the first data card in each data set. 

Any number of data sets may be stacked behind the initial data control cards. 
The first data set of each job setup must use: the OPR, CAL, and BASE cards; 

the T1 and AT cards; and, the SKYTEMP cards. A /* card may be used to trun- 
cate succeeding data sets in any one of three places: after the OPR card; 

after the BASE card; or, after the AT card. 

"Radome Is On" Data Card 


The "Radome Is On" card (example follows) contains the Log Id., the valid 
data code (VDC), the disk address, and the encoded antenna elevation angle of 
the disk data set that contains the antenna loss values needed when computing 
the radome loss. 


60 

o 

-j 


o 

£ 


« 

a 

4) 

U 

3 

a 


at a> 

0 r-4 

fi 00 

S5 

•3 c 

o 

V tI 
ai 4J 
•O CO 

o > 
u a> 
C rH 

w w 


r Tasumr is 'an" nm oo 


iDiJlfisE 






- 4 


t ' •*■■•** v ' ♦' 

l ODOOOIOflO 01> 00 0 6 Oil’ . tltOOOflQ.OO 00. , 

: ( M 1 1 1: m nnNMjJimil 

t V n 11 1 U V » i ij.^i ii i j ill y j 1 1 u 1 1 h i w 1 1 i*ij i *i 

2Z2I72222 22222222212222222 222222 21222 22222 2 2 2 2 2 2 2 2 2 2 2 2 22 2 2 2 2 2 ^22 2«2 2 2 22 24^.^ 
3 3 3 3 3 3 3 3 3 3 3 3 1 3 2 3 3 3 3 3 *3 3 J 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 ) 3 3 3 3 il 3 3 #4.3 3 3 3 33333 31 312#. : 
444. 4. 44444*4444 « 44*<v044444*4444444444‘4444 4W4 4 4l44 4*4V44^*^474 44i4tl^W^ 

* 5 5 55 5 5 5 5 55 5 5 5 5 5 5 5 535 555 5 55 55 5 55 5 5 5 5 5 55 5 55 5 55 5 55SS.5S5S55S 55$i555566955555ltjU^ 

* 66 05 6666 6 6 6 6 6 6 6 6 6 6 6 6666666 6 66 6( 66 66 6 6 6 6 6 60 6 ( 6 6 6 6 6 6 6 6 6 6(6 6 6 6.6 I MMt$6£^£l* , 

^ , I? t h ? p W? r i 1 ? 7 m 7 7 7 *1 rnTiii fit ifi Iniilli 11 111111 1 ii ) ? Jl) )fl 7 imfi 1 

?mm. er MM itata# ww L 

" J 6 55393 5 p MSI6l9i9 999 9 9 S9l95S99#l5 1 9 9 91 9 9 9 9 9 1 9 9 9 It 9ff 9 9 

1 1 i 1 1 j > 1 * i«^ 8 n » an»naa»aana»ana«j4»aaaa^yc«)«<wa oa^a> 


HIM] 


"RADOME IS J " CARD - OUTPUT BY MFMR 3 "RADOME IS OFF" JOB 

The Disk Addresr on the "Radome Is On" card is used to access engineering 
data, and it is therefore not possible to av^age several antenna losses (as 
in PMIS 3) through manipulation of the data deck. 


page M3. 9 


134 


REFERENCES 

Blume, H-J. C. and C. T. Swift, n S-Band Radiometer Remote Sensing 

of Sea Surface Temperatures," Proc. 1972 URSI/USNC Conference, 
Williamsburg, Va., Dec. 12-15, 1972 (USNC-URSI, Nat'l. Acad. 
Sciences, 2101 Constitution Ave. , N. W. , Washington, D. C. 
20418) . 

Carver, K. R. , "Remote Sensing Using Microwave Radiometry , " 

Proc. 1973 IEEE Southeast-Con , Louisville, Ky., April 30 - 
May 2, 1973 (IEEE, 345 E. 47th St., New York, N. Y. 10017). 

Carver, K. R. , Wm. Cooper and J. F. Paris, "Local Variations in 
Radiometric Sky Temperature Using Radiosonde Data," 

Proc- 1975 URSI/USNC Conference, Urbana, 111., June 2-5, 1975. 

Fogarty, Wm. G. , "Total Atmospheric Absorption at 22.2 GHz." 

IEEE Trans. Ant, and Prop. , AP-23, May, 1975, pp. 441-444. 

Hach, J.-P., "A very Sensitive Airborne Microwave Radiometer 

Using Two Reference Temperatures," IEE Trans. Micro. Theory 
and Techniques , MTT-16, September, 1968, pp. 629-636. 

Hansen, R. C. , Microwave Scanning Antennas, Vol, I , Academic 
Press (New York), 1964, pp. 24-46. 

Kerns, D. M. , "Correction of Near-Field Antenna Measurements 
Made With an Arbitrary But Known Measuring Antenna," 
Electronics Letters , Vol. 6, No. 11, May 28, 1970. 

Kraus, J. D., Radio Astronomy , McGraw-Hill Book Co. (New York), 
1966, p. 261. 


PBECEDING PAGE BLANK NOT FILMED 



135 


REFERENCES (Continued) 

Kraus, J. D. and K. R. Carver, Electromagnetics , 2nd ed. , McGraw- 
Hill Book Co. (New York), 1974, pp. 621-622. 

Mentzer, C. A., L. Peters, and R. C. Rudduck, IEEE Trans ■ Antennas 
and Prop. , AP. 23, March, 1975, pp. 153-159. 

Paris, J. F., "A Program for Computing the Brightness Temperature 
of a Clear Atmosphere from Radiosonde Data," Tech. Rept. 
LEC/HASD No. 6490.21.068, Lockheed Electronics Co., Inc., 
Aerospace Systems Div. , Houston, 1971. 

Paris, J. F., "Sky Brightness Temperature Error Analysis for 

Microwave Sensor Calibration," Job Order 75-415, Lockheed 
Electronics Co., Inc., Aerospace Systems Division, Houston, 
Feb. , 1975. 

Sullivan, W. T., "Variations in Frequency and Intensity of 1.35 - 
centimeter H2O Emission Profiles in Galatic HII Regions," 
Astrophys. J. , vol. 116, 1971, pp. 321-332. 

Van Vleck, J. H. and V. H. Weisskopf, "On the Shape of Collision - 
Broadened Lines," Rev. Mod. Phys. , vol. 17, 1945, pp. 227-236. 

Yu, J. S., R. C. Rudduck and L. Peters, IEEE Trans. Antennas and 
Prop. , AP-14 , March, 1966, pp. 138-149. 


